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CHAPTER 1 
INTRODUCTION 

A. Cosmic Microwave Background Background 

The existence of a blackbody background radiation as a permanent 
remnant of the flash of radiation of the "big bang" was first suggested by 
Gamow and his collaborators (Gamow 1948, Alpher and Herman 1958). Using 
Friedman's (1922, 1924) solution to the field equations of General Relativity 
for an expanding isotropic homogeneous universe, Gamow and his associates 
attempted to explain the present element abundances m terms of nuclear 
reactions occurring during the first few seconds or minutes of the expansion. 

They showed the following (Alpher and Herman 1950, Gamow 
1949, 1953, Alpher, Herman, and Gamow 1967): 

1. At early epochs radiation would be in thermal equilibrium and 
thus would have a thermal spectrum. More important, at that time it would be 
the dominant component of the total energy density. 

2. As the universe expanded the radiant energy decreased as the 

inverse fourth power of the scale factor. On the other hand, the matter rest energy 

density decreased only as the inverse third power of the scale factor ; hence 

eventually it became the dominant component of the total energy density. 

4 

3. At an epoch when the temperature^/as approximately 10 °K , the 
radiation decoupled from the matter. In spite of the complexities of this de- 
coupling process the distortion of the blackbody spectrum of the radiation by 
this process would be small. Peebles (1968) has since treated this process 

in some detail and has shown this to be essentially a consequence of the fact that 




the energy density of the radiation at the time of recombination of the primeval 
9 

plasma (~10 ev/baryon) is very much greater than the 10 ev/baryon released 
by hydrogen recombination. 

4. Following the decoupling, the radiation would appear to an ob- 
server who was a rest with respect to local matter as Doppler shifted blackbody 
radiation. The well known result that Doppler shifted blackbody radiation be- 
comes blackbody radiation with a lower temperature implies that the radiation 

would still have a thermal spectrum. 

3 4 2 

5. Synthesis of He , He , and H is possible during the first few 

minutes of the expansion of the universe It is important to note here that 

Gamow's original aim was to account for the genesis of the majority of elements, 

4 

but there now' seems little doubt that elements heavier than He cannot be so 
produced (Fermi and Turkevich 1950, Wagoner, Fowler and Hoyle 1967, 

Peebles '1966). They are presumably made m s<,ars (Burbidge, Burbidge, Fowler, 
and Hoyle 1957). 

Gamow realized that the radiation would persist to the present epoch, 
but unfortunately he did not suggest any attempt to detect it experimentally. It 
was not until the work of Dicke and his collaborators (1966) who were looking 
for observable consequences of an expanding universe that the background rad- 
iation was reconsidered,* and they specifically constructed a specialized radio- 

* Independently Doroshkevich and Novikov (1964)exammed existing observational 
data in an attempt to determine if there was any evidence for the existance of the 
radiation considered by Gamow. Unfortunately, due to a misreading of a paper 
by Ohm (1961) which described the results obtained at the Bell Telephone Labor- 
atories with an absolutely calibrated corner horn reflector and a quiet receiver, 



telescope for its detection. However, the radiation was iirst detected by 
Penzias and Wilson (1965), when it was pointed out to them by Dicke that the 
excess antenna temperature of unknown origin of their exceptionally quiet 
microwave radiometer might be due to the radiation which the Princeton group 
was seeking. The second observation of the radiation followed within a year, 
and was that of Dicke's co-workers Roll, and Wilkinson ( 1966) Their measure- 
ments as well as others made subsequently, are shown m Figure 1. The best 
present estimates suggest ~2. 7°K for the temperature of this radiation. 

One of the important predictions of the theory is that the radiation 
should have a spectrum corresponding to that of a blackbody. For a temper- 
ature ~3°K this would peak at about a millimeter. m wavelength. As can be 
seen from Figure 1, the thermal character of this radiation has been confirmed 
over more than two decades m frequency in the long wavelength portion of the 
spectrum However, since most of the energy density lies at the short wave- 
length portion where measurements have not yet been made, it is of great 
importance that observations be extended to .this spectral region If the 
thermal character of this radiation is. eventually confirmed, it will provide 
very compelling, and perhaps conclusive evidence in favor of the expanding 
universe 

they came to the conclusion that there was no evidence for its existence. Penzias 
and Wilson (1965) eventually used the same instrument to first detect the rad- 
iation. 
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B. Observational Difficulties 

Direct observations have not been made at short wavelengths because 
of the ‘presence of atmospheric H O and O absorption lines. These increase 
in both strength and number as one goes to shorter wavelengths. Their resulting 
opacity causes the atmosphere to radiate with an intensity which is large compared 
With that of the background. In addition, the blackbody background intensity 
begins to fall off exponentially with increasing frequency beyond about a milli- 
meter in wavelength. Thus it appears that beyond \ «3mm, for all practical 
purposes, the "radio window" is closed, and one must place radiometers above 
the atmosphere m order to make direct background measurements. 

Fortunately, it is possible to- circumvent the problem of atmospheric 
opacity by using observations of optical spectra which are due to absorption by 
the interstellar molecules CN, CH, CH+. An analysis of the absorption spectra 
of these molecules can be used to determine to what extent the molecules , are 
being excited by the background radiation. A measure of this excitation can 
then be used to infer the intensity of the existing background. 

This dissertation describes how this is done. The resulting measure- 
ments made by the use of these molecules are shown in Figure 1 > where they are 
labeled CN, CH, CH+. 

C. Interstellar Molecules 

The optical interstellar absorption lines due to the interstellar mole- 
cules CN, CH, CH+, were discovered by Adams and Dunham (1937, Dunham 
1941, Adams 1941, 1943,1949). Following the suggestion of Swings and Rosen- 
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feld (1937) that the line at X 4300. 3A might be due to CH, McKellar (1940) 
succeeded in identifying this feature with CH as well as the line at X 3874. 6 A 
with CN. At that time he predicted the presence of several other lines of CH 
as well as the presence of the R(l) line of CN at X3874. OA. Acting on the sugges- 
tion of McKellar i Adams succeeded in observing the additional lines of CH as 
well as the faint CN R(l) line. Douglas and Herzberg (1941) then produced 
CH+ in the laboratory, and positively identified the features at >,4232. 6A and 
X 3957. OA as due to that molecule. 

It is surprising to consider that McKellar' s (1941) estimate of 2. 3°K 
for the rotational temperature of the interstellar CN m front of C Ophiuchi was 
probably the first measurement of the cosmic microwave background — twenty- 
eight years ago.' His temperature measurement was made from Adams' visual 
estimates of the absorption line strengths, so his result was crude. At the time, 
though, he attributed little significance to this result stating (McKellar, 1940) 

". . . the 'effective' or rotational temperature of inter- 
stellar space must be extremely low if, indeed, the con- 
cept of such a temperature in a region with so low a den- 
sity of both matter and radiation has any meaning. . . " 

His 2. 3°K measurement is also mentioned in Herzberg's (1959, p. 496) Spectra 
of Diatomic Molecules , suggesting that it ", . . has of course only a very re- 
strictive meaning. " 

The starting point of our work was a crucial suggestion by N. J. Woolf 
based in turn on a discussion by McKellar (1940) which preceeded Adams' dis- 
covery (McKellar 1941) of CN X 3874. OA. He suggested that the absence of 
excited state CN lines placed a severe limit on the temperature of background 
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radiation As we shall see, not only do these molecules set upper limits to this 
radiation, but m all probability they may also be used to effect an actual measure- 
ment of it 

In addition to our work on this problem (Thaddeus and Clauser 1966, 
Clauser and Thaddeus 1969, and Bortolot, Clauser, and Thaddeus 1969) similar 
work has also been pursued by Field and Hitchcock (1966), while Shklovsky (1966) 
made an early suggestion that McKellar's (1941) observation was a consequence 
of the background radiation. 

D. Purpose of D issertatio n 

The central purpose of this dissertation is to extract -the largest 
feasible amount of information on the short wavelength spectrum of the micro- 
wave background radiation from a number of spectra of interstellar molecules 
that were available m the summer of 1966 The work may be divided into roughly 
two parts, an observational part and a theoretical part In the observational part, 
we first show how one calculates the background intensity from the spectra of 
the interstellar molecules (Chapters 2 and 3). Second, in order to utilize the 
existing spectrograms, we develop new techniques of spectrophotometry (Chapter 
4, Appendix 4) Third, we present the results of the application of these tech - 
niques to the problem at hand, and the resulting upper limits to and measurements 
of the intensity of the background radiation (Chapter 5). 

On the theoretical side, we first consider what assumptions concerning 
the location of the molecules will be necessary, such that (1) the molecules will 
set upper limits to the background intensity, and (2) the molecules will yield 
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reliable measurements of the background intensify. It will be seen that the use 
of the molecules to set intensity upper limits places only minor restrictions on 
the molecular environment, however, their use to make intensity measurements 
requires the absence of alternate molecular excitation mechanisms (Chapter 2). 

We therefore provide an analysis which considers possible excitation 

i 

schemes that are consistent with our present understanding of the conditions found 
in the interstellar medium We find that a necessary assumption for the mole- 
cules to yield reliable intensity measurements is that they reside in a normal 
H I region (Chapters 6,7, and 8) 

After a review of the existing evidence that the molecules do 'reside 
in an HI region, we proceed to present new evidence to further substantiate 
this contention. This wilL consist of (1) the observed invariance of the excitation 
of the interstellar molecules, and (2) the results of new 21cm H I observations 
in the direction of these molecular clouds (Chapter 9). 
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CHAPTER 2 

THE USE OF INTERSTELLAR MOLECULES TO OBTAIN UPPER LIMITS AND 
MEASUREMENTS OF THE BACKGROUND INTENSITY 


A Excitation Temperature of Interstellar Molecules 

The method which we exploit to obtain information concerning the 
short wavelength spectrum of the backgiound radiation is the familiar one used 
to find the temperature of molecules and their environment from molecular 
spectra (Herzberg 1959). It is a general tool widely used in astronomy to 
determine the temperature of planetary and stellar atmospheres, and finds 
considerable application in the laboratory to the study of flames, rocket 
exhausts, etc . 


/• 
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f 

\ 


\ 




0 or B 
STAR 



HI CLOUD SPECTROGRAPH 


A typical observational situation is illustrated schematically above 
The intensity in absorption of lines in a molecular electronic band is proportional 
both to the square of the transition matrix element and the population of the lower 
level of the transition As long as the molecules are m thermal equilibrium 
(or at least their rotational degrees of freedom are in thermal equilibrium) the 
relative populations of the rotational levels are in turn given by the usual Boltz- 


mann expression 
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n. g. 

— = — exp [ -(E. - E.) / kT 1 (2.1) 

n. g. x 3 

where 

n. is the population of level i and 
g is the statistical weight of that level 
E is the energy of that level 

It is now clear that the relative intensity of optical absorption lines in the band 
is a function of T alone , and may be used to determine T if the rotational term 
scheme and relative transition matrix elements of the molecule are known, as is 
usually the case. Even if thermal equilibrium does not hold, it should be obvious 
that equation (2. 1) can be used to define an "excitation 1 ' or "rotational" temp- 
erature T.. with respect to any two levels i and i. 

13 

Consider now possible causes of this excitation; it could be due to 

either radiation or collisions. If radiation which has a thermal spectrum with 

temperature is the only (or dominant) cause of the molecular excitation, 

then on purely thermodynamic grounds T.. = T . In this case only the emission 

lj B 

and absorption of photons with frequencies 

V. =(E. - E.) / h = / 2tt (2.2) 

i] i 3 13 

m the molecular frame of rest are involved in the excitation; and only the portion 

i 

of the radiation spectrum over the narrow Doppler width 

_g 

A v. / v « 10 
i] 1J 

of the transition is involved in the excitation of the molecules. Thus if the 
radiation does not have a thermal spectrum, but over these intervals it has the 
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same intensity as that of a blackbody temperature T , then the relation 

13 

T = T will hold: 

B ij 

On the other hand if the excitation is collisional, or is due to a 
mixture of radiative and collisional processes, T will generally set an 
upper limit to the intensity of radiation at Thus we see that a measurement 

of T will yield useful information on the background radiation density. 

B . Astronomical Situation 

In the astronomical situation of interest to us the molecules are 
located in a tenuous interstellar cloud, while the source of radiation against 
which absorption is being observed is an O or B star at a typical distance of 
100 pc. Since the interstellar lines are very sharp, and for the molecules, at 
least, always weak, the largest telescopes and Coude spectrographs are required 
to detect the lines. In fact, virtually all the molecular observations have been 
made with three instruments; the Mt. Wilson 100-inch, the Mt. Palomar 
200-inch, and the Lick Observatory 120-inch reflectors. 

There are several reasons why interstellar lines have so far only 
been investigated against early stars. Perhaps the most important is the high 
intrinsic luminosity of these objects, which allows high resolution spectroscopy 
to be done in reasonable exposure times to distances of up to several hundred 
parsecs. For much smaller path lengths, the probability of intercepting an 
interstellar cloud which possesses appreciable molecular absorption is very 
slight. There is the additional factor that O and B stars have been so recently 
formed that they usually lie near the galactic plane, and thus m the vicinity (or 
at least the direction) of interstellar matter Also an important technical con- 
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sideration is that the spectra ol these stars are comparatively leatureless, many 
ol the spectral lines are due to ionized atoms and the lines are usually wide 

* i 

due to the high temperature and/or rotation of the star. In marked contrast to 
the situation with cool stars, there is thorelorc little likelihood that the inter- 
stellar lines will be blended or obscured bv stellar features. Purely -from the 
technical pome of view wo thus see that early stars are ideal objects to use as 
light sources for long range optical absorption spectroscopy of interstellar matter. 


C. Calculation of T . 

ij 

It is straight forward to calculate T from- the observed absorption 

intensities of lines in a molecular band The observed optical depth t , of 

lr 

an absorption line, which is due to a transition from level 1 1 to level 1 '. , is 
given by 

t . - 8nv , NS , /. (Uhcg ) (2 3) 

ii* n' i i i 1 i 

where S , is the line strength This is defined as the square of th’e electric- 
dipole matrix element summed over the degenerate magnetic states and polanza- 
tions of the transition (Condon and Shortley 1951, p. 98). 


S = X) i 'imjji Cj I jmJ ~ (2.4) 

J m m q 4 J 

i I 

where M is the electnc-dipolo momenl of the molecule and C =,/4tt /3 

lq 

^lq ’ w ^ iere i ^ c are the normalized spherical harmonics This 

summation may be expressed in terms of the reduced matrix element (Edmonds 


1968, equation 5. 4. 7). 
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S =^ 2 1 ( i H c (l) |1 j >| 2 (2.5) 

ij ~ 

where C (1)=/4tt/3 Y(l) is a spherical tensor operator of order one. 

Thus for two optical absorption, lines, of nearly the. same wavelength, 
that are due to transitions to levels i' and j' from low lying levels i and j 
respectively, 


r 


Here we have defined 



= s exp 



( 2 . 6 ) 


s 


n' 



( 2 . 7 ) 


The calculation of this quantity for the specific molecular systems of interest 
will be taken up in Chapter 3. Equation (2. 6) may be rewritten to give the 
excitation or rotational temperature 

hv.. 

T ij = kln(s7r) (2 ' 8 

It should be noted in passing that since only the logarithm of the 
optical depth appears m equation (2. 8), the rotational temperature may be well 
determined even though the optical depths are poorly known. 

If the interstellar lines are optically thin, the optical depths are 
proportional to the equivalent widths of these lines (area of absorption line for 
unit continuum height). Fortunately, most of the interstellar lines considered 
in this study are weak so that the correction for saturation is small, and good 
approximate results may be obtained by taking r equal to the equivalent width 
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ratio. However, accurate temperature measurements will require a correction 
for saturation via a curve of- growth for the clouds which contain these molecules. 

A description of these corrections appears in Chapter 5. 

D. Why Molecules ? 

We must now consider which kinds of atoms or molecules found m 
the interstellar medium are useful as radiant thermometers, for measurement 
of the background radiation. What is required to make such measurements 
possible is an atom or molecule which has (1) a low lying level with energy 

separation from the ground state kT « 2cm \ for a blackbody radiation 

» 

temperature T ?«3 0 K] (2) an allowed (electric-dipole) transition connecting 

this level with a lower level which has an observable population and (3) allowed 

optical absorption transitions originating m each of these levels. Clearly these 

optical transitions must be mutually resolvable, and of course, ‘to be'observable 

from the ground they must have wavelengths X > 3000A. The required term 

diagram for such a species is shown in Figure 2. 

We can see that there is nothing m princsple which restricts us to 

use molecules. However, atomic fine- structure separations are typically much 

greater than 2cm . For example, the lowest lying fine-structure level of any 

interstellar atom so far observed from the ground belongs to Ti II and this is 

97cm 1 above the ground state. On the other hand, hyperfine levels have 

separations which are too small, and are connected by magnetic -dipole transitions 

-1 

- the H I hyperfine separation is only . 032cm . 

The rotational and/or fine structure energy level separations within 
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the electronic ground states of the three molecular species observed optically in 

the interstellar medium - the CN, CH, and CH+ free radicals - have energies 

not very much greater than 1 88cm 1 . The J = 0 - l energy level separation 
2 

of the X S (v=0) electronic ground state of CN is especially favorable with a 
separation of only 3. 78cm \ Partial term diagrams for these molecules are 
shown in Figure 3 Precise molecular data is listed in Table 1. 

2 + o + 

The R(0) absorption lines of the (0, 0) band of the B 2 - X 2 

electronic system of CN at X3874. 0A and X3874. 8 A arrise from the J = 1 

level. All of the CN lines are easily separated with a good Coud£ spectrograph 

It should be noted, however, that even in the best spectrographs the lines are not 

fully resolved, l. e. the observed line width , ~0. 02 a, is not the true Doppler 

line width of the interstellar cloud, but is determined instead by the resolution 

of the instrument, and is typically » 0 05 A . 

The level separations of the CH and CH+ molecules are sufficiently 

great that no molecules are observed to be in an excited level. However, these 

two molecules are still of great value to the problem of the background radiation, 

since the absence of the excised states lines can be used to set upper limits to 

the radiation intensity at very short wavelengths. 

The lowest excited state in CH+ is the J = 1 level, 27. 9cm 1 above 

the J = 0 ground state, and the corresponding R(0) and R(l) transitions of the 
11 

(0, 0) band of the X S - A r electronic system occur at X4232. 5A and X4229. 3A 

2 

The CH electronic ground state configuration is a F_ state, inter- 
mediate between Hund's cases (a) and (b). Accordingly, transitions are allowed 
from the KJ = 1 l/2 ground state to both the fine-structure level KJ = 1 3/2, 
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which lies at 17. 8em 1 above the ground state, and to the K J = 2 3/2 excited 

-1 

rotational level 66. 8cm above the ground state. The strongest optical trans- 
ition from these three levels are the R (1), and R (1), and the R (2) lines in 

£ _L /L 

2 2 o o 

the (0,0) ‘band of the X If- A A electronic system at X4390. 3A, X4303 9 A 

o 

and 1,4296. 6 A respectively. A detailed treatment of the relative intensities of 
the optical transitions of CH appears in Chapter 3. 

The J = 2 level of the CN electronic ground state is similarly useful 
for setting an upper limit to the background intensity Since it is 7. 6cm ^ above 
the J = 1 level, its populations will be quite small. For this reason, the R(2; 
line originating m this level has not yet been detected, but it is probably only 
about a factor of three below detection (Bortolot, Clauser, and Thaddeus 1969) 
and a further extension of the techniques developed in this dissertation might 
yield its detection. 

On the other hand, without an enormous extension of present ob- 
servational techniques, the CH R (L) and the CH+ R(l) lines are below the 

a 

threshold of detectability if their temperature ~ 3°K. 

E. Stat is tical' Equilibrium of an Assembly of Multi-Level Systems Interacting 
with Radiation 

When all its components are considered, departure from thermal 

equilibrium is one of the most conspicuous features of the interstellar medium 

4 

It contains, for example, dilute 10 °K starlight, H I atoms with a kinetic temp- 

4 

erature of ~ 100° K, 10 °K H H protons and electrons, grains which are thought 
to be in this vicinity of 20°K, and high energy cosmic rays. It is obvious that 
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in 'ts totality.it cannot be characterized by a unique temperature. 

At first glance there would, therefore, appear to be no reason to 
assume that interstellar CN should single out the microwave background with 
which to be in thermal equilibrium. However, for the rotational degrees of 
freedom of CN it will be shown that this is probably the case. 

In order to discuss the statistical equilibrium of interstellar molecular 
systems, it is probably wise to be specific about a few conventions and defin- 
itions that will continually appear during the discussion. 

Thermal Equilibriu m: An assembly of systems is in thermal equilibrium when 
all of its parts are characterized by the same temperature. 

Statistical Equilibrium: If an assembly of systems with energy levels is in 
statistical equilibrium when the populations of these levels have time-independent 
populations 

Due to the long time scale of interstellar processes (and the short 
time scale for observations) there are often grounds for assuming statistical 
equilibrium to apply, even though extreme departures from thermal equilibrium 
are known to exist. All of the assemblies which we will be considering 1 . e. 
clouds of interstellar molecules, will be m statistical equilibrium. This may 
be seen from the fact that an isolated polar molecule such as CN, CH, or CH+, 
interacting background radiation, will achieve statistical equilibrium with this 

r* 

radiation typically in less than 10°sec . On the other hand, it will exist m 
the interstellar medium for at least 10 ^ sec, before it is destroyed by radiative 
dissociation. Thus the assumption of a constant number of molecules m statis- 
tical equilibrium is clearly a very good one. 
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Excitation Temperature: Boltzmann's law .states that if an assembly of systems 

with energy levels E. is in thermal equilibrium, then the populations in these 

levels will be distributed according to Boltzmann's law - equation {2. 1) - 

which may be used to define T in terms of the level populations. If there is 

thermal equilibrium at temperature T,then all T = T. Alternately if thermal 

equilibrium is not established, equation (2, 1) may be used to define T as the 

i] 

excitation temperature of the levels i, j. * 

Brightness Temperature : Boltzmann's law, when applied to the radiation field, 
yields Planck's law for the density of unpolarized isotropic thermal radiation 


u(v) = 


o , - 3 -3 
8rrhv c 

/hv 

exp fe 


-3 -1 

(erg cm Hz ) 


(2. 9) 


and u„ = u( v_). If the spectrum of the radiation is not thermal, it is conven- 
tional to use equation ^2. 9) to define the brightness temperature T (v..). Since 

B ij 

u^ is a monatomc function of T^(v_), the latter is frequently used m this 

case as a measure of the radiation density at frequency v . 

Principle of Detailed Balan ce: The principle of detailed balance states that the 

transition rates for an elementary process and Its inverse process are equal 

This Is just due to the invariance of transition matrix elements under time - 

reversal (see for example Tolman 1938, p. 521)*. Hence if the transition rate 

for a given process R from level i to level j is given by R. , then 

i] 


* Although the principle of detailed balance does not strictly hold in a relativistic 
theory including spin (Heitler 1954, p. 412) - for our purposes it may be consid- 
ered correct. 
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g.JR.. = g R.. exp-j/-(E. - E.)/kTl (2.10) 

J Ji ii] i J 

where T is the kinetic temperature of the exciting agent. In the following 

analysis we will make the important assumption that kT E.-E., which is true 
for the molecular levels of interest if the kinetic temperature of the colliding 
particles is 125° K or higher. Then we may set the above exponential equal to 
one. A temperature m the vicinity of 125°K is the usual one derived for i. 
normal H I region(Dieter and Goss 1966). The suggestion that shielding from 
the interstellar radiation field by grains may yield significantly lower temper- 
atures in dense, dark H I regions does not apply to the clouds in which CN, 

CH, and CH + have been observed, since these clouds are reddened less than 
1 magnitude. 


The statistical weight factors arrise in (2 10) from the fact that 
the elementary process in question occurs from one state to another, but the 
labels i and j refer to levels Following Condon and Shortley (1951) we refer 
to a level i as the set of g^ magnetic states with the same angular momentum 
and energy E. 

Einstein Coefficient A.. (sec 1 ): A is the probability per unit time that a 

ij i] 

system will spontaneously emit radiation of frequency v . and make a transition 

O 

from level i to level j. Hence the rate of such transitions for n. molecules in 

3 i 


level lis'n.A... 

i i] 

3 

Einstein Coefficient B.i sec erg cm Hz): B u. is the probability per unit 

13 1} 

time that a system will absorb radiation from a field of energy density u„ at 

frequency v and make a transition from level j to level l. For n. molecules 
ij J 

in level j the absorption rate is then just n B..u. . For consistency with the 

3 i] i] 

principle of detailed balance B u. must then be the rate of stimulation of emission 

ij H 

of radiation at frequency v .. by the radiation field u... For n molecules in 

ij i] i 
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level i the rate for this process is equal to n.B..u... For a two level system 

1 

in equilibrium with thermal radiation then 


if we wish to identify 


n , B. u 

_i = p - __JUJ — 
n. s ii A..+B .u 
1 i] 13 13 


T_ = T.. then we must have 
B 13 


B . 
3 31 



and 



„ . 3 
Srrhv.. 
11 


3 

c 


( 2 . 11 ) 


( 2 . 12 ) 


(2. 13) 


F. Equations of Statistical Equilibrium: " 

' Let us now consider an assembly of multilevel systems such as 

diatomic molecules in an interstellar clouq. We shall assume the molecules to 

be excited by background radiation and other processes as well. These might 

be collisions, fluorescent cycles through higher electronic state, etc. , but we 

shall make no initial restrictions on them, at least at first, specifying them only 

by R , as the rate of transitions induced by them from level J to level J T . 

For statistical equilibrium of the n molecules in level J, we 

J 

require 



0 =-n 


E 

J'=0 


R 


JJ 1 


+ E 

J’=0 




(2. 14) 
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J’<J 


- Tf C n j (A jj> + b jj t V ■ n j |B j'J u JJ’ ^ 

J ,— o 

+ £ t n jT (A J’J + 

J' >J 


(2. 14 
cont'd) 


For interstellar molecules in the presence of non-thermal radiation 

and additional interactions T. . will then be a complicated function of the various 

13 

transition rates given by the solution of equations (2. 14). Fortunately, we do 
not need to solve these equations in general, but can obtain results for three 
specific cases relating the molecular excitation temperature to the radiation 
density at the frequency of the microwave transition. The first two will apply 
to CN and CH+ , while the third will apply to CH. 

In the treatment of cases I and II, which refer to diatomic molecules 
with 2 electronic ground states, the levels will be denoted by their total 
(excluding nuclear) angular momentum J. In the treatment of case III which 
refers to a molecule such as CH with both spin and orbital angular momentum 
in its electronic ground state, we will, for simplicity of notation, refer to the 
levels KJ = 1 l/2, 1 3/2, and 2 3/2, as 0, 1, and 2, respectively. 


G. Case I : 

Let us first consider the simplest case of interest, that of a cloud 
of interstellar diatomic molecules which posses a 2 electronic ground state 
(so that the allowed rotational transitions satisfy the simple selection rule 
AJ = ±1) interacting only with background. Assuming that higher multipole 
transitions may be neglected, we will show that ^ j ~ T^(Vj t ^ j). 
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It may be seen in this ease, that different pairs of adjacent levels 

may have different excitation temperatures, and that these temperatures are 

equal to the brightness temperatures at the frequencies of the corresponding 

transitions. The various T will not, however, be equal when the spectrum 

ij 

of the radiation is non-thermal. 

The proof is straight forward. If the allowed rotational transitions 

satisfy the selection rule A J = ± 1 and all R . = 0, then the linear homogeneous 

i] 

equations (2. 14) reduce to 


dn 

~dT 


+ n j+l^ A J+ 1, J +B J+1, J U J+1, j* " D J B J, J+1 U J+1, J 


(2 15) 


This system of equations may be solved' by Inspection by noticing that the first 
and second lines both vanish identically if one takes 


n j+l _ ^ B J, J+1 U J+1, J 

n -r %+l, J A Ti . +B Ti , T u__ 

J J+1,J J+1,J J+1,J 


(2 16) 


This is just a generalization of the solution (2. 11) for a two level system Thus 
upon combining equations (2. 9), (2. 11), (2. 12) and (2. 13), we have that 


T = T (v ) 

J+l, J E> J+l, J' 


(2- 17) 


H. Case II: 

Consider now the addition of other interactions (e. g collisions) to the 
situation of Case I. For the problem of immediate interest, we know from ob- 
servation that the J = ,0 and J = 1 rotational levels are the most populated levels. 
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and that these two populations are not inverted (or for the case of CH+ , only 

the J = 0 level is observably populated). 

We will show that and T„_ will in general set an upper limit 
—10 21 

to the brightness temperatures T (v ) and T (V ). Hence if the CN R(2) 

B 1U Jo 

and the CH+ R(l) absorption lines are not detected (assuming that the corres- 
ponding CN R( 1) and CH+ R(0) lines are observed) an upper limit to their optical 

l 

depths will set upper limits to the radiation densities at the frequencies of the 
CN (1 -» 2) and the CH+ (0 - 1) rotational transitions. 

The proof for this case is somewhat more involved. If the allowed 
rotational transitions satisfy the selection rule A J - ± 1, but R_ - 0 (we make 
no assumptions concerning the selection rules for the additional processes), then 
equations (2. 14) specialize to 


dm 




® = ° = Sjj ' + jS nj ' Rjj ' 


" n J (A J, J-1 U J, J-l* + n J-l B J-l, J U J, J-l 


(2. 18) 


+n j+l (A J+l, J +B J+1, J U J+1, n J B J, J+1 U J+1, J 

Using equation (2. 10) we can rewrite equations (2. 18) for J = 0 


f H 1 W 

— + e + j'=i 0J ' \ W 

n .JQ ' 

0 A 10 B 10 U 10 


(2. 19) 


and for J = 1 
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n 


n. 


n_ 


Vi + 


A 10 +B 10 U 10 “ E 01 U 01 + J f^ 1 H ,J 1 


cr n 

& 1 J y 

g J ,n i / (2. 20) 


A + B u 
21 zl 21 


Consider first equation (2. 20). 'We are given' (from observation) that n /n < 

J 1 

g /g for all Jn 2. Hence the individual terms of the summation are all 

ci X 

positive so that nVn - and T - will be an upper limit to T (v_. A ) 

Also note from equation (2. 14) that we can redefine u^q such that R = 

= 0. Let us denote this by u and in the same respect define 

^1 0 a? 10 . 'We note that - 




n. 


* 


10 


( 2 . 21 ) 


since we have introduced no negative terms into (2. 19) by this .proceedure. 

We are now in a position to consider equation (2 20). Introducing 


U 10 and ■= io * nto ^is e< 4 ua fci 0il we have 


B u’ 

!g _ . + Jii° V^LO " 

n =21 


n \ » 

i -/Ks » 


1/ J=2 


1J 


V_J 

g A, 


A 21 + B 21 U 21 


By our previous reasoning and by the use of inequality (2. 21) we find that 

n /n ^ E , thus T n , will be ah upper limit to T ( v o-.) 

2 1 21 21 B 


I. Case III: 

Finally we proceed to consider the case of an interstellar cloud of 
CH molecules excited by both radiation and collisions. Again we appeal to 
observation and notice that the lowest level is by far the most populated. 
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CH is a diatomic molecule with both spin and orbital angular momen- 
tum in its ground state. Its coupling scheme is intermediate between Hand’s case 
(a) and (b), hence the selection rules A J = 0, ±i, A K = 0 ± 1, ±2, apply. 

Since the spontaneous decay rates increase with the cube of the 

♦ 

energy for higher levels, molecules in these levels will very rapidly decay back 
to the ground state. Also since the populations of these higher levels are very 
small pur calculation need only consider the three lowest lying rotational-fme 
structure levels. These are shown in Figure 3a. There are the levels KJ = 

1 l/2, 1 3/2, and 2 3/2, which for brevity we will label 0, 1, and 2, respectively. 

We will show that in contrast with the case of molecules with only 
AJ = ± 1 transitions permitted, the fact that the transitions have a non-zero 
branching ratio prevents us fropa always using the population ratio to set an 
upper limit to the radiant intensity. We will find that n^/ and n^/n^ will 
not simultaneously set upper limits to the radiation density at frequencies 
and \\ , but that (1 + &,) n,/n n and (1+ ej n/n_ will set upper limits, 1 . e. 

v oi 1 1 U Z Z U 


? 10* (1+e l ) T* and e 20“ 1+ V^- 


& and e are dependent upon the conditions of the molecules 1 location, and for 

1 2r 


- 5 ~12 

e fn 10 and e « 10 
1 & 


conditions typical of the interstellar medium 
Thus for all practical purposes, the level ratios themselves may be used to 
define the radiant intensity upper limits. 


Since all these levels are connected by radiative transitions, following 

the proceedure used in Case II we will absorb all of the Ik. into u*. > u„, 

and m a similar fashion define f 1 sf; .It will be noted that the f (or 

13 13 "13 ** 



KJ Labe! 



Figure 3 a . Transitions A try to keep levels 0 and 1 in thermal equilibrium. 
Transitions B and C depopulate level 1 and populate levels 2 and 0. The 
population of level 1 may become depressed if the rate of transitions B 
becomes very large with respect to that of transitions A. 
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are not necessarily related except in the case of thermal equilibrium of the 
whole system - then ? 1 o^2i = W 

Specializing equations (2. 14)' to this case, we can write the express- 
ion analagous to (2. 19) and (2. 20).. - 


n_ 


n. 


" S' 


10 


r 


10 


and 




20 


t< ti 
5 21 s 10 


- 1 


D, 


(2. 23) 


n 2 

— - V 
n Q { 20 


V, 


20 


V' 


1 - 


20 


? 2r io 


Dr 


(2. 24) 


with the denominators given by 


D 1 = 


B u' 

01 10 

B u’ 

02 20 


V, 


20 


? 2r io 


B u' 
12 21 


B u' 

10 + _W10_ + x 


B 12 d '21 


(2. 25) 


D _ B 02 U, 2Q 
2 B_u'. 


1 + 


*10 \ B 20 U '20 

4- + 


r 


20 


B 12 U 2l/ B 21 U '21 * 21 ^ 10 


(2. 26) 


01 10 

Unfortunately expressions (2,23) and (2. 24) cannot be simultaneously 
positive, i. e. either is bound above by n^/n Q or 5^0 by n 2^ n 0‘ Tllus 
in marked contrast to the cases of CN and CH+, it is possible to depress the 
population of the first excited CH level by additional interactions. This can be 
seen from Figure 3a. If we arbitrarily increase the rate of 1^2 transitions 


while keeping 0 -> 1 and 0 -*• 2 rates constant, we can effectively remove mol- 
ecules from level 1 and transfer them to level 2, from which a portion will decay 
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to level 0. This process will effectively decrease ‘n^/n^. It will occur if the 

1 2 rate becomes large with respect to radiative 1 0 decay rate. 

Let us consider equations (2. 23 - 2. 26) slightly further. If 

A / (B u' ) is large, corresponding to a small rate B ’ u* , then the 
XU X<« tjj. ±.£ 

denominators D and D will be large and n /n and n / n will set very 
X ^ 1 U U 

good upper limits." This tepm is always large unless we have an enormous flux 
at . For CH ^ « 200p the strongest suggested source of radiation 

at this wavelength is thermal radiation from interstellar grams (Partridge and 
Peebles 1967). For the expected brightness of these, 

A 10 = lA l - 11,5 < 2 - 27 > 


so is bounded by n^/n^ (1 + 10 ). 

The same term also appears in the denominator D multiplied by the 

z 

factor B n 1 / (B u’ ) . Since u’ is now bounded, all that is required is a 
Uo iZv ux XU t XU 

minimum energy density af v , e. g. that radiated by grains, and D will 

^u - z 

also be large. Taking the maximum energy density at allowed by observation, 
and the minimum at provided by grains, then 



Hence we also are safe in using n /n (1+10 .) to bound ^ . 

£ U ciXj 

Nor can we call o.i collisions to significantly weaken these upper 


limits. In order to weaken the upper limits, we require a selective excitation 


rate for R with no associated increase in R_ and R , otherwise, the 

- 1 UX U Z • 

populations of at least one of the levels, 1 or 2, would have an observable 
population. 
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However, excitation of molecular rotation by collisions (of the first 
land) with projectiles whose kinetic energy is very much greater than the, 0-1 
and 1-2 level spacing is not a selective process. Thus, collisions cannot pro- 
vide the selective 1 -> 2 excitation necessary to depress the K J - 1 3/2 level 
population. In addition, even with the largest excitation cross sections due to 

the longest range forces (Coulomb), charged particle densities of at least 100 
-3 

cm would be required for there to be any effect at all As we will see later, 
there is specific evidence 'that the charged particle densities are nowhere near 
this high. 

The only way to significantly weaken our upper limits is to have 
narrow band radiation at ~200ja which selectively excites the KJ = 1 3/2 - 
2 3/2 transition of CH. The temperature calculated for this radiation must 
be at least as high as 12° K. Although such "line radiation" cannot be totally 
ruled out by present observational work, it seems highly unlikely that such 
radiation might exist m the interstellar medium with an intensity this great. 

It must be sufficiently narrow band that it not excite the CH XlSOp. transition 
and the CH+ X 359 p transition. 

Thus it appears that we are quite safe m general in assuming that 
n /n and n /n set safe upper bounds to §' and and that the temper- 

I U 6 U JLU 2jU 

atures corresponding to these population ratios provide good upper limits to the 
brightness temperatures at and v 0 ». 

JLU 

J. Sufficient Conditions for Thermal Equilibrium to Hold 

From equations (2. 14) and (2. 15) it can be seen that for molecules 
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with a 2 ground state, n./n^ will yield -an actual -measurement <of T^v^) 

if all R « B u . where is the intensity necessary to produce 2 7°K 
ij 01 10 10 

excitation. At present we need only concern ourselves with CN, since this is the 
only molecule m which an excited state population has been observed. * 

The spontaneous decay rate may be written (Condon and Shortley 
1951, p. 98) in terms- of the line strength s j ^ and rotation constant B as 


J,J-1 


„„ 4 3 
647T v 

3hc 3 (2J+l) 


S 


4 3 
512tt B 


J,J-1 ' 3h 


2J+1 J,J-1 


(2 29) 


Application of Edmonds’ (196 8) equation (5. 4. 5) to equation (2. 5) yields 


l j) =/j 


(2. 30) 


m terms of the Wigner 3j symbol, and hence 


A 


J,J-1 


4 3 
5 127T B 

3h 


J 4 2 
M 


2J+1 


(2. 31) 


From a study of the optical stark effect of CN, Thomson and Dalby 
(1968) have recently determined the CN dipole moment to be 1. 45 ± . 08D. 

This agrees well with the predictions of molecular orbital calculations (Huo 1967) 
which indicate that the dipole moment is expected to be in the range 1. 2 - 2. 5D; 
and more or less confirms the work of Arpigney (1964) who deduced (J. = 1.1D 
from an analysis of cometary spectra. 

The rotation constant B Q of CN has been determined to high precision 

by an analysis of the violet band to be B^ = 1. 8909cm (Poletto and Rigutti 1965) 

* Expressions are derived in Appendix 2 which may be used for evaluating the 
corresponding Einstein coeffcienN Tor CH. 
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Using these values of and p, equation (2. 33) yields for CN 

= 1. 1889 x 10 5 sec 1 
A = 11.413 x 10“ 5 sec" 1 

ai 

What we must show m order to claim that the CN is in thermal 
equilibrium with the background radiation is that the R„ are much smaller 
than 


B oi u io 


f] 

Sr 


10 


fhv 


= 8 69 x 10 


exp | 


10 


kT 


- 1 


(2. 32) 


where u^ is the intensity al necessary to produce an excitation temperature 
T = 2. 7°K. 

For most quantum mechanical processes the excitation rates fall off 
very rapidly with increasing AJ since the dipole term in the interaction poten- 
tial is usually the dominant one, and has selection rules AJ = 1 Thus it will 
be sufficient for our purposes to simply show that for a given process, R 
« B 01 u 1() . We will see in Chapters 6, 7, and 8, that this is probably so if 
the CN is in an HI (non- ionized) region of space of reasonable density, but 
that it is not necessarily so if the CN resides in an H II (ionized) region. 
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CHAPTER 3 

OPTICAL TRANSITION STRENGTH RATIO 


A. CN and CH+ 

In the last Chapter we saw that the rotational temperature may be 

t 

calculated in terms of r, the observed optical depth ratio, and s , th'e theoretical 
transition strength ratio defined by equation (2. 8). 

For lines in the same band s is independent of the electronic os- 
cillator strength and Franck- Condon factor, and often depends upon only the 
angular momentum quantum numbers of the participating levels, K, K’, J, and 

J ! . The CN violet system is a good example of this. Although this transition 
2 2 

is Y) Z) > the p- doubling produced by the electron spin is too small to 
be resolved in interstellar spectra, and on the grounds of spectroscopic 
stability (Condon and Shortley 1951, p. 20) we are justified in calculating 
intensities (and assigning quantum numbers m Figure 3 as our use of J for K 
implies) as though the transition were ^Z/ - * ^Z) • 

Thus when a molecule is well represented in either Hund's case 
(a) or (b), as are both CN and CH+ , the S , are proportional to the Honl- 
London factors (Herzberg 1959, p. 208), and we have 


J. + J„ + 1 

i i* 


J. + J. t + 1 


(3. 1) 


However, for the case of CH which is intermediate between these 
two representations, the calculation is somewhat more involved. It requires 
a transformation of the matrix elements of the electric-dipole moment operator. 
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calculated in the pure Hund's case (b) representation, to the "mixed" represen- 
tation in which the molecular Hamiltonian matrix is diagonal. 

To effect this transformation, we need the matrices which diagonalize 
the molecular Hamiltonian matrix; hence, we also calculate this Hamiltonian 
matrix in the Hund's case (b) representation, and then calculate the transformation 
matrices which diagonalize it. As a check we notice that the resulting line strength 
ratio agrees with equation (3. 1) when we specialize our results to the case of a 
rigid rotor. 

B. Matrix Elements of the Molecular Hamiltonian 

Since the rotational energy of CH is large with respect to its fine 
structure, even for the lowest rotational levels, we calculate the molecular 
Hamiltonian in Hund's case (b). In this representation the appropriate coupling 
scheme is (Herzberg 1959, p. 221) 

K = N + A , J=K + S (3.2) 

where in units of h 

A is the component of electronic angular momentum along the 

A-' 

internuclear axis, 

N is the molecular rotational angular momentum, and 
S is the spin angular momentum of the electron. 

Appealing again to spectroscopic stability we ignore X. - doublets and hyperfme 
structure since these interactions, like p- doublets, are unresolved m inter- 
stellar optical spectra. 

Following Van Vleck (1951), the molecular Hamiltonian is 
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K = B (K 2 A 2 ) + AA S , (3. 3) 

~ z 

where A is the spin-orbit coupling constant, "S , is the spin angular momentum 

z 

2 . 

projected along the internuclear axis, and B = ti /2I , where I is the molecular 

moment of inertia. Ther term in S , mixes together states of different K when 

z 

the fine structure is appreciable with respect to the rotational energy and effects 
the transition away from Hund's case (b) towards case (a). 

The matrix elements of this Hamiltonian for the Hund's case (b) 
coupling scheme are calculated in Appendix 1 to be* 

(AK'SJ | V | AKS-J) = B[K(K+1) -A 2 ]6 k , k 

+ AA(-l) A+K+K+S+J r(2K'+l)(2K+l)l 1//2 ( j S K'j /K' 1 k\ 

Jl K S j \A 0 h) 

If the 3j and 6j symbol are evaluated algebraically, this single 
expression yields both Van Vleck's (1951) equation (2 4) for the diagonal, and 
his equation (25) for the off-diagonal matrix elements, but with the opposite 

i 1 

sign convention for the off-diagonal terms. ** 

* Van Vleck (1951) has calculated these matrix elements by exploiting the an- 
omalous commutation relations of the angular momentum components referred 
to the molecular frame. He showed that the problem corresponds to the familiar 
case of spin-orbit coupling m atoms However, since his results yield the 
opposite sign convention from ours, and are of little value for calculating intensities, 
we have recalculated them here by more standard means. 

** A misprint occurs in Van Vleck's equation (25). An exponent of 1/2 has been 

2 2 l/2 

omitted from the last factor, which should read [(K + 1 ) - A ] . The correct 
expression appears in the earlier paper of Hill and Van Vleck (1928). 
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Here selection rules may be seen implicit in the properties of the 
3j and 6j symbols. The selection rule A K = 0, ±1 for the molecular Hamil- 
tonian in case (b), for example, follows from the requirement that the rwo 
K', 1, K of the 3] symbol m equation (3. 4) satisfy the triangular rule. 

For a molecule intermediate between Hund's case (a) and (b), in 
general, K is no longer a good quantum number m the representation where the 
Hamiltonian of equation (3. 4) is diagonal. We may, however, label a given state 
with the value of K which it assumes when the molecular fine structure constant 
A vanishes. We will distinguish K by a carat when it is used as a label m 

A 

this way The stationary state is then written as \ A KSJ ) , and the unitary 
matrix which transforms from case (b) to the diagonal frame as U~ . This 
unitary matrix will be used to transform the dipole moment operator into the 
representation m which the molecular Hamiltonian is diagonal, and may be 
calculated from (3. 4) by a standard matrix diagonalization proceedure. 

For the diagonalized molecular Hamiltonian matrix itself, we have 

(AKSjj y j AKSJ > - X) U^, K u| k <AK'SJ| 3 f j AKSJ) (3.5) 

K'K 

Inserting into this formula the coupling constants for CH of Table 2, and 
diagonalizing the matrices for the upper and lower electronic states of CH, 
we have the partial energy level scheme shown m Figure 3. 

C. Line Strengths 

As before, we take the line strengths to be proportional to the 
absolute square of the electric dipole moment operator reduced matrix element 
S =p 2 j (AKSJa’ 1 C (1) |1 AK'SJ'co')! 2 , (3.6) 
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where a represents all other non-diagonal quantum numbers. 

This reduced matrix element has been evaluated m Appendix 2 for 
the CH coupling scheme 


< AKSJal' C (1) I 1 A'K'SJ'a') 


& D k. K . u kk 


, ^.A+S+J'-fl T -1/2 

(-1) f (2K’+l)(2K+l)(2J'+l)(2J-i-l) y /A 


If J 8 ljS<«l C lq |a.(> 


K 1 K'\ ,2 
A q A7 1 


(3. 7) 


The usual symmetric-top selection rules for electric-dipole transitions , 

AK = 0, ±1, A J = 0, ±1, are implicit in the 3] and 6j symbols (but note that 
the transition to case (a) allows AK = ± 2 transitions). 

When both the upper and lower states are diagonal m the Hund's case 


(b) representation, and we let the electronic matrix elements ( a 1 | Y [a) 

* 

equal one, expression (3. 7) reduces to that of the Honl-London factors (Herzberg 

» l f 

1959, p. 208). 

When equation (3 7) is specialized to the case of a rigid rotor (CN 
and CH+) with K' = J', K = J, and S = 0, and with the electronic matrix element 
again equal to one, using Edmonds' (1968) equation (6. 3. 2) we have 


I <AJ|j C(l) !' A\P>f = (2J r +l)(2J+l) / J l j'\ ( 

\-A q Ay 

where q = A- A' and \ q ( ^1 

With the molecular constants found m Table 2, the resulting values of s, cal- 
culated from equations (3. 7) and.(3. 8)* for the coupling schemes of CN, CH, 


* 


in agreement with (3. 1) for CH+ and CN 
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and CH+, are found in Table 1. 
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CHAPTER 4* 

TECHNIQUES OF SPECTROPHOTOMETRY AND PLATE SYNTHESIS 

A. Available Plates and Spectrophotometry 

The optical lines of interstellar CN are characteristically sharp and 
weak, and occur in the spectra of only a handful of stars. The strongest line, 
R(0) at X3874. 6 A , has an equivalent width which is seldom as strong as 10mA, 
while the neighboring R(l) and R(2) lines are weaker by a factor of ~4 and 
~100 for an excitation temperature ~3°K, 

The state of the art in high resolution astronomical spectroscopy for 
the detection of weak lines has not greatly improved over the past thirty years - 
the weakest line which is detectable with a reasonable degree of certainty on a 
single spectrogram is ~ 1 - 3 mi Hence it was clear at the outset of this work 
that a new technique would have to be developed if we were to attempt to detect 
the CN R(2) line, as well as to measure accurately the R(l) line. 

At the beginning of this work the best available spectragrams of 
interstellar molecular absorption lines were those of Adams and Dunham from 
Mt Wilson, Miinch from Mt. Palomar, and Herbig from the Lick 120-mch 
telescope. Since long exposures are the rule for these spectrograms, and 
observing time on large telescopes with good Coude spectrographs is somewhat 
difficult to obtain, it appeared that the most straight forward way to improve the 
signal- to-noise ratio was to add the many available individual plates, and to filter 
out the high frequency gram noise from the composite spectra. 

Direct superposition of the spectrograms on the densitometer 
("plate stacking") is not practical when many spectra are to be added together, 
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or when the spectra are dissimilar. In addition, this proceedure does not allow 
one to independently weight the individual spectra. -The more flexible method 
of computer deduction of the data'was felt to be best.- This technique involved 
digitization of the spectra, and their calibration and synthesis on a digital 
computer. * The spectra could then be filtered of gram noise numerically in 
the manner suggested by Fellgett (1953) and Westphal (1965). The whole system 
turned out to be very flexible so that the large heterogeneous group of plates 
could be added without difficulty 

Most of the plates were obtained from the extensive collection of 
high-dispersion spectrograms stored in. the Mt. Wilson Observatory plate files. 
There were obtained primarily by Adams and Dunham with the Coud6 spectrograph 
of the 100-meh telescope and its 114" camera, and were used by Adams (1949) 
in his classic study of interstellar lines. The spectra of the two ninth mag- 
nitude stars BD+66°16?4 and BD+66°1675 were taken by Munch (1964) with 
the 200-inch telescope. The synthesis for C, Ophiuchi -includes six spectra 
obtained by Herbig (1968) with the Lick 120-inch telescope using the 160" camera, 
and were used by him in his recent study of the interstellar cloud m front of 
Q Ophiuchi. (The latter were also analyzed by Field and Hitchcock [19641 in 
their discussion of the microwave background. ) We are indebted to Drs. Munch 
and Herbig for the loan of these plates. 

* Similar techniques have been used. by Herbig (1966) for a small number of 
similar plates 
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Data such as the emulsion, exposure time, etc. , for the individual 
plates were copied for the most part from the observing notebook, and are 
listed in Table 3. Also in this Table are the relative weights (determined 
from statistical analysis of grain noise as described below), which the indi- 
vidual spectra are given in the final synthesis. 

It is evident from the Table and from the weights that they by no 
means represent a homogeneous collection. The weights were determined by 
a visual estimate of the gram noise r. m. s. level on tracings which were nor- 
malized to the same signal level. The normalized tracings were then combined 
with the usual prescription of weighting them inversely as the square of the 
noise level. 

These plates were scanned at. the California Institute of Technology 
on a Sinclair Smith microdensitometer, the output of which was fed across the 
Caltech campus via shielded cable to an analog- to- digital converter. This was 
connected to an IBM 7010 computer which wrote the digitized signal on magnetic 
tape, and subsequent analysis of the data was .done on an IBM 360/75 computer 
at the Goddard Institute for Space Studies (GISS) in New York City. 

Figure 4 shows the setup of equipment at Caltech The output of 
the densitometer was digitized at a constant rate of 100 samples per second. 

The lead screw was not digitized but also ran at a constant speed. The projected 
acceptance slit width and lead screw speed were adjusted to suit the dispersion 
of the plate and the width of the narrowest spectral features being studied, 
which were about 0. 05 A. Thus for the Mt. Wilson 114" camera plates, 

o 

(dispersion = 2. 9A/mm) which comprise the bulk of the data, 
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the plate advanced at 5. 6mm/min (16. 3A/min). For the narrowest lines this 
is about five line widths per second. In other words at least 20 data points 
were obtained per resolution element, which is ample for our purposes. 

- The projected acceptance sht width was adjusted to the dispersion 
of the plate being scanned to approximately 1/6 of the full line width of the 
interstellar lines. Thus for plates from the 114" camera the projected slit 
width was about 2. 9 p, which leads to negligible signal distortion. 

Four or more densitometer traces were usually taken m the vicinity 
of the three wavelength regions of interest, X 3874A X42321, and X4300A . 

Two were required for the comparison spectra on either side of the stellar 
spectrum, at least one was required for the stellar spectrum itself, while a 
single scan perpendicular to the dispersion of the spectrograph was necessary 
to record intensity calibration wedge bars. 

In order to establish a reference wavelength common to the stellar 
and comparison spectra, a razor blade was placed across the emulsion as 
shown in Figure 5. A scan started on the razor blade and moved off the edge 
onto the spectrum as shown in the Figure, with the edge producing an initial 
step in the recorded spectrum. In addition to providing a wavelength fiducial 
mark, the step provided a measure of the response time of the electronics and 
the optical resolution of the system, and in the case of the intensity wedge, 
provided a dark-plate level. 

B. Synthesis Programs: 

The computer program for plate synthesis was divided into three 
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sub-programs, each subprogram was processed separately with the. input for 
a given sub-program being 'the output.tapes from the previous one. ' There were 
named Sys 3, Sys 4, and Sys 5. 

’Sys 3 read the 7040 tapes and data such as the transverse lead screw 
position, which had been punched manually onto cards, and scanned this data for 
errors. Preliminary plots of each vector versus index were made after it had 
been adjusted to start at the fiducial mark step (caused by the razor edge), so 
that the index of a given digitized point m the vector was proportional to its 
corresponding distance from the razor edge. Output of plates corresponding 
to different molecules were sorted onto different tapes so that processing of a 
tape after this point was for a single spectral region The basic purpose of 
Sys 3 was to pre-process as much information as possible to implement the 
smooth running of Sys 4. 

Sys 4 took the left adjusted vectors for -each plate and converted 
these to two (or more if there were more than one input stellar spectrum vector) 
parallel vectors of intensity and wavelength (shifted to the molecule's rest frame). 
This was done in two parts using an IBM 2250 CJRT display console m conjunction 
with the GISS computer. The first part used the wedge vector to convert the 
recorded plate transmissions to intensity. The calibration density wedge bars 
had been put on the plate at calibrated light levels, the values of which were 
stored in the computer. 

The wedge bars from the density wedge vector appear as rectangular 
peaks as shown in Figure 6. The limits of these rectangular bars were deter- 
mined m Sys 4 by displaying the wedge vector on the ‘2250 console, and picking 
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put the ends of the rectangular -peaks with .a light pen. The computer then 
averaged the transmission over each wedge bar. Since the two sides of the 
W.edge w,ere put on the plate at different intensities,- the wedge values on one 
side were scaled up or down to interleave them between adjacent points on the 
other side. The composite curve was then fit with five terms of a Laurent series, 
a typical result of which is shown m Figure 7. The series was finally applied to 
each point m each spectral vector converting the plate transmission to intensity. 

The second part of Sys 4 calibrated the index scale in terms of wave- 
length Each comparison spectrum was displayed on the 2250 Display Console 
and the centers of selected lines were located with the console light pen. The 
corresponding wavelength of each line was then plotted on the 2250 screen as a 
function of index. Since the dispersion is highly linear over the ~10A traced 
in the vicinity of each interstellar line, these points appeared on the screen m 
a straight line. A straight line least-square fit was also displayed on the screen, 
which aided m the correction of any line misidentifieation The output of Sys 4 
was finally two parallel vectors of intensity and corresponding wavelength appro- 
priate to the laboratory frame. 

In general, high resolution Coude spectrographs are rather well 
matched to the emulsions used in astrophysical spectroscopy, and grain noise 
which is "high frequency" with respect to sharp interstellar lines is not con- 
spicuous, In order to extract the maximum amount of statistical information 
from the data at hand, however, a considerable amount of thought was given to 
the problem of numerically filtering the final spectra. A description of the 
statistical analysis appears in Appendix 4. The final result of the analysis 
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was to prescribe a’filter function to be used to filter the calculated synthesis. 

Thus, Sys 5 was designed to take the successive outputs from. Sys 4 
and to add these together with the appropriate statistical weights. It then filtered 
out the ''high frequency" noise by convolving the synthesis with the above filter 
function. The final output of Sys 5 for ea'ch star and molecule was three vectors: 

1. Wavelength in the molecular rest frame 

2. Unfil’tered spectrum 

3. Filtered spectrum 
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CHAPTER 5 
RESULTS OF SYNTHESIS 

A, Curve of Growth Analysis 

Figures 9, 10, and 11, show the results of adding together', in the manner 
just described the available spectra of interstellar CN, CH, and CH+ for a 
number of stars. As just mentioned the wavelength scales iri these plots are 
those of the molecular rest frame, and the final spectra except for £ Persei 
CN and £ Ophiuchi CN have been numerically filtered. The number of in- 
dividual spectra included in each synthesis and the corresponding vertical mag- 
nification factor are indicated on the left side of Figure 9. 

Table 4 lists the visual magnitudes, MK spectral types galactic 
latitudes as well as the equivalent widths W and optical depths t , for the inter- 
stellar lines measured. The uncertainties of W m all instances result from a 
statistical analysis of the gram noise, and represent a level of confidence of 
95. 5% (see Appendix 4). Also all upper limits set on the strengths of unobser- 
ved lines represent the more conservative confidence level 99. 7%. 

For £ Ophiuchi the equivalent widths of the interstellar molecular 
lines listed in Table 4 are systematically smaller by ~ 30% than the values 

obtained by Herbig (1968a) from plates taken at the 120-inch Coude We note, 
however, that Dunham (1941) reported W - 6mA for the CN it(0) line and 
W = 14mA for the CH R (1) line in £ Ophiuchi, and O. C. Wilson (1948) reported 

A 

W = 16 ± 2mA for the CH+ R(0) line in £ Ophiuchi, both from Mt. Wilson plates 
and both in good agreement with out values of 6 62, 13. 4, and 19. 6 mi re- 


spectively. 
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Since most of the spectra from which our values are derived were 
also obtained with the 100-inch telescope about thirty years ago, this gives us 
some confidence that the discrepancy is not a result of an error in calibrating 
or adding together the spectra, but represents instead a systematic difference 
between equivalent widths obtained with the 100-mch Coude then and the 120-inch 
Coude now. * Our main interest, however, is only in the ratio of line strengths, 
which m the limit of weak.lines is independent of systematic uncertainties of this 
kind, and it is accordingly unnecessary to explore this point further. But we 
emphasize that the uncertainties listed in Table 4 can only be taken seriously m 
a relative sense, and are liable to be considerably less than absolute uncertainties. 

t 

Also included m this Table for completeness are several late B and 
early A type stars. whose rotational temperatures have been estimated by other 
workers. 

The calculation of the optical depths in Table 4 requires some know- 
ledge or assumption concerning the true shape and width of the interstellar molec- 
ular lines. In general the lines we are considering are not quite resolved, and 
it is thus necessary to make several assumptions concerning the curve of growth 
Fortunately the CN lines are usually so weak that the correction for saturation 
is small, and our final results are not very sensitive to these assumptions. 

* It is interesting to note, however, that Adams (1949) was led to suggest that 
the molecular absorption spectra may be time variant. 
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We assume a gaussian line shape in all cases (see Appendix 4) and 
use the curve of growth due to Stromgren (1948) and Spitzer (1948), and the 
tabulation of Ladenburg* (1930). For % Ophiuchi we adopt for the linewidth 
parameter b = 0. 61km/sec, which is the value b = 0. 85km/sec found by Herbig 
(1968a), for interstellar CH+ in the spectrum of this star, corrected for the 
apparently systematic difference m the measured equivalent widths of the CN R(0) 
line. In lieu of any other information we assign b = 0. 85km/sec tothe CN 
lines of £Persei, since this star lies at a distance comparable to £ Ophiuchi, 
and is also well off the galactic plane. ** (We have also found with the 300-foot 
transit radio telescope at Green Bank, that the 21cm profile in the direction of 
this star closely resembles that m the direction of ([Ophiuchi - see Chapter 9. ) 

The signal strengths for the CN lines of the remaining stars shown 
m Figure 9 scarcely justify correction for line saturation: we simply adopt 

b = lkm/sec for all remaining stars except BD+66°1674 and BD+66°1675, to 
which we assign b — 5km/sec, on the basis of their great distance and low 
galactic latitude. 

The excitation temperatures of interstellar CN may now be calculated 
from the ratios of the optical depths listed m Table 4 from equation (3. 4). 

* This was rechecked and a finer tabulation made. 

** Further work should be done to obtain a good curve of growth for this star 
due to its importance here. 
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B. CN (J - 0 -» 1) Rotational Temperature and Upper Limit to the Background 
Radiation at X 2. 64mm 

Table 4 now shows the results of the application of equation (3. 4) to 

11 -1 

this data, taking v Q = 1. 134 x 10 sec ( X = 2. 64mm) corresponding to the 
j = 0 -* 1 transition in CN. The listed value of T Q (CN) results from the ratio 
of the optical depths listed in Table 4. The entry labeled AT lowers T^(CN) 
by an amount which is at most comparable to the uncertainty due to the plate grain 
noise and , therefore , our final results as mentioned above do not critically 
depend on the assumptions that were made concerning the curve of growth 

lb get an upper limit on the radiation intensity we need consider only 
£ Ophiuchi , where the molecular absorption lines have the greatest signal-to-noise 
ratio. For this star the weighted mean of the values found from the line strength 
ratios JR(1)/R(0) and P(1)/R(0), given respective weights of 4 to 1, the square 
of the corresponding values of s (see Chapter 3), . yields 1^(2. 64mm) ^T^CN) = 

2. 74 ± 0. 22°K. In terms of the intensity of radiation this becomes 

— X5 -2 — 1 — X -1 

1^(2. 64mm) ^3. 42 x 10 erg cm sec ster Hz 

C. Upper Limit to Background Radiation at X = 1. 32mm - 

The observed absence in the spectrum of £ Ophiuchi of the CN R(2) 
line at X 3873. 369A places an upper limit to the background intensity at 
X = 1. 32mm. From the spectrum of Figure 9 it is found that the equivalent 
width of the CN R(2) line is less than 0. 44mA, to a level of confidence of 99. 7%. 
The curve of growth analysis described above then yields t< 0. 031, and the . 
rotational temperature found from equation (3. 4) and the R(2)/R(l) ratio is 
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T (1. 32mm) <T (CN) < 5. 49°K 
B 4&1 


and 


-14 -2 -I -1 -1 

1^(1. 32mm) <2. 75 x 10 erg cm sec ster Hz 

which is still a factor of 8. 3 above the intensity of 2. 74°K thermal radiation, 
but a factor of 1. 6 below the intensity of 2. 74°K greybody*radiation. This 
upper limit is shown m Figure 12. 


D. Radiation Upper Limit at X = 0. 359mm 

The observed absence of the CH+ R(l) line at X 4229. 341 A in the 
spectrum of Q Ophmchi synthesized from 19 spectra of this star and shown in 
Figure 11, allows an upper limit to be placed on the intensity of radiation at 
the wavelength of the CH+ J = 0 -> 1 transition (0, 359mm). 

As with CN the optical depths listed in Table 4 were calculated using 
Herbig's (1968a) value b = 0. 85km/sec for the linewidth parameter, adjusted 
to take into account the apparently systematic difference between our and his 
equivalent widths. 

For the CH+ R(1)/R(0) ratio it was found that s = 3/2, = 8. 35 x 

lO^^Hz and 

T_ (0. 359mm) < T (CH+) <9. 01 K, 

B , 


or 


-14 -2 -I -1 -1 

1^(0. 359mm) < 10. 1 x 10 erg cm sec ster Hz 


* i. e. the extrapolation of the Rayleigh- Jeans intensity I = 2kTv 2 / c 2 to short 
wavelengths. Radiation from an optically thin medium such as dust can approximate 
a spectrum of this type. 
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13 + 

B. Detection of Interstellar C H 

13 + 

Also shown in Figure 11 is the faint trace of the R(0) line of C H 
originally sought by Wilson (1948) and again by Herbig (1968b). A statistical 

o 

analysis yields for the equivalent width 0. 6 ±. 3mA to a confidence level of 95%. 

12 + 

Together with the R{0) line of C H this yields for the isotopic abundance 
ratio along the line of sight C^/C^ - 62 m agreement with the value 

— CiCt 

+55 

82 __ independently found by Bortolot and Thaddeus (1969) and also with the 
-15 

terrestial value 89, but much greater than the equilibrium value for the 

CN bi-cycle (Caughlm, 1965). This observation strengthens Herbig's (1968b) 

earlier conclusion that the material m the Q Ophiuchi cloud did not come from 
13 

C rich stars, but raises the question of how one explains the apparent agree- 
ment of this value with the terrestial one. 

Bortolot and Thaddeus (1969) suggest that the coincidence of the 
13 

terrestial C abundance with that of the interstellar abundance is due to the 

13 

origin of the solar system out of an interstellar medium with a fixed C abun- 

* 13 

dance. This notion is consistent with the coincidence of the atmospheric C 

abundance of Venus with the terrestial abundance ( Connes et al 1968); however it 

is inconsistent with the recent observations of Lambert and Maha (1968) who 

13 12 

find the solar abundance ratio to be C /C = 150 ±30, a value which exceeds 
the terrestial ratio by nearly a factor of two. * 

Lambert and Malia believe that their observation represents at least 


.13 


* There is considerable dispute on the question of the solar C abundance 
(Cameron, 1969). 



56 


a definite upper limit to the solar C abundance. Since there is no obvious 

» 

13 

process which might increase the solar C abundance, further work confirming 

their observations will represent strong evidence favoring the suggestion of 

Fowler et al (1961) that the planetary abundances are a product of spallation, 

13 

and is otherwise unrelated to the interstellar C abundance. 

F. Radiation Upper limits at X= 0. 560 and X= 0. 150mm 


Figure 10 shows the result of the synthesis of seventeen spectra of 

C Ophiuchi m the vicinity of X = 4300 A, the strongest lines of CH. The observed 

absences of the CH R (1) and R 2 (2) lines at X= 4304. 946A and 4296. 636A allow 

upper limits to be placed on background intensity at the wavelengths of the 

CH KJ = 1 1/2 - 1 3/2 (X - 0. 560mm) and KJ = 1 l/2 - 2 3/2 (X= 0. 150mm) 

transitions. We again use b = 0. 85mm (corrected). 

For the CH KJ = 1 l/2 *■» 1 3/2 transition s was calculated m Chapter 

11 

3 to be 1. 525. Taking ^ = 5. 37 x 10 Hz ( X = 0. 560mm), equation 

(3. 4) and the data of Table 4 then yield 

Tg(0. 560mm) ^ ^ x ^ (CH) < 6. 33°K 


-14 _2 -1 ~i 

I^jO. 560mm) ^3. 91 x 10 erg cm sec ster Hz 

Finally for the CH KJ = 2 3/2-*l l/2 transition s = 0. 820 and 
12 

V 2 3/2 1 1/2 = 2 ' 00 x 10 Hz = 150mm). Hence 


T b ( 0. 150mm) * T 2 3/2j 1i/2 (OH)< 28. 5°K 


1(0. 150mm) < 4. 2 x 10 ^erg cm % sec ^ ster ^ Hz . ^ 




57 


G. Discussion of Upper Limits 

At 150/Li the most important source of background might be the 

contribution from young galaxies ( Low and Tucker 1968, Partridge and 

5 

Peebles 1967). However the CH upper limit at 15 O^i is~10 greater than the 
expected integrated brightness of these objects and as a result is of little value. 

The constraints which the other upper limits impose on the background 
radiation are shown m Figure 12. Except for the point at 2. 64mm, they fall well 
above the blackbody curve at T = 2. 7 °K, and all are considerably below the 
corresponding greybody curve. They immediately represent much of the infor- 
mation we possess on the radiative content of the universe at these wavelengths. 

It bears final emphasis that as upper limits, all of these points are virtually 
independent of any considerations concerning the location or environment of the 
molecules. 

Shivananden, et. al. (1968), have recently reported the detection of 

a background flux in the wavelength interval from . 4 to 1. 3mm with a rocket- 

-9 ~2 -1 

bo rne wide-band detector. They detected a flux of 5 x 10 W cm ster 

-14 

which yields a mean intensity over this frequency range of I = 9. 6 x 10 
-2 -1 -1 

erg cm ster Hz . This flux is represented in Figure 12 by the area 
under the line AA'. Since its mean value is nearly three times greater 
than the CH upper limit which falls right zn the center of their frequency range, 
it may be concluded that either their observation is of a background with highly 
irregular spectrum (eg. spectral lines), or it is of a flux of quite local origin 
which does not extend as far as the major interstellar cloud m front of Q 


Ophiuchi. 
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The important question which is unanswered at this point in the 
analysis is what is the meaning of the observed J ~ 1 level population of CN. 
In the remaining Chapters we will try to show that, in all probability, this 
population provides a good measure of the radiation intensity at X 2 64mm, 
and the intensity at this wavelength is equal to that of a blackbody temperature 

t b “ T 10 (CN) - 
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CHAPTER 6 

ALTERNATE RADIATIVE -EXCITATION MECHANISMS 

In addition to direct rotational excitation of CN by microwave photons, 
radiative excitation is possible to both higher vibrational and electronic states 
by infrared and optical photons. The characteristic lifetimes of these 
higher states are very short, and virtually all of the molecules will be found m 
the ground electronic and vibrational states. However, indirect rotational 
excitation is possible by means of electronic and vibrational fluorescence. A 
CN molecule can be excited to a high vibrational or electronic level by an optical 
photon (from dilute starlight for example) and re-emit a photon while making a 
transition back to a rotational level different from the one m which it started. 
Hence, this cycle can effectively populate higher rotational levels, and contri- 
bute to the observed rotational temperature. Typical cycles are shown m Figure 
13 . 

Such indirect processes are potentially very serious alternate mech- 
anisms for the CN rotational excitation, in that they might provide an explanation 
for the observed rotational temperature invariance in eleven stars. This observed 
invariance provides the strongest single piece of evidence that the CN is not 
excited by processes other than the cosmic microwave background. (See Chapter 

9 ) 

The explanation for the invariance in terms of fluorescence could 
come in either of two ways. First, if the excitation is produced by galactic 
starlight, the inherent uniformity of this radiation might account for the observed 
invariance. Second, the process of CN formation m the interstellar medium 
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is not well understood (see, for example, Dipter and Goss 1968), and it seems 
plausible to suggest that the CN is only formed at- certain distances from a star, 

, and thus would always be subject to the same radiation field. In this respect* 
Munch (1964) has proposed that the CN is being flashed off interstellar grains 
by the advancing ionization front at an H I - H n boundary. 

However, it will be shown in this Chapter that, (1) vibrational 
fluorescent excitation is negligible by more than six orders of magnitude. (2) 
Electronic fluorescence is insufficient to account for the observed CN rotational 
excitation by at least three orders of magnitude for galactic starlight. (3) If 
the CN is close enough to a star, perhaps the star used to observe it,for fluores- 
cence to contribute materially to the rotational excitation of the ground state then 

* v 1 "* 

the photo-dissociation rate will be so high that the molecules will be destroyed at 
a rate much faster than any suggested mechanism for their production. In add- 
ition, the Idea that the molecules are cireumstellar is inconsistent with 
out present ideas concerning the location of the molecules. 

A. Criterion for Eluorescence to be Negligible 

We saw In Chapter 2 that a sufficient criterion for a process to be 

negligible is for ft „ „ <<A , where J" and J"+l represent 

J } J "r J. o J 

rotational levels in the electronic ground state ** e”. 

* The suggestion by Merrill (1942, 1946) that the CN might be cireumstellar was 

dispelled by Adams (1949, see Chapter 7). 

\ 

** Unless otherwise specified, e n will be used to denote either an electronic 
or a vibrational state. 
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Fluorescent excitation involves first a transition from rotational 

level J" m the ground electronic and vibrational state e", to a level J' in a 

higher electronic or vibrational state e'. This will occur at a rate given by 

u „ ,B „ , Tl . Subsequently, the molecule decays to a rotational level 

e , J , s e tl J l e'J", e'J' 

J in the electronic ground state e". Since the excited vibrational or electronic 
state lifetimes are very short, in general, 

A e'J',e"J >> U e'J , ,e ,, J nB e"J",e , J r - (6 ' 1) 

Thus for a given band system, the total rotational excitation rate is given approx- 
imately by 



a 


.u 


B 


e'J', e"J n e'J, e' J" e"J",e T J' 


( 6 . 2 ) 


where a is the branching ratio.of e’J -» e"J electronic decays 

6 J j 6 J 

(including cascades through an intermediate electronic state) to all other possible 
electronic decays, a will generally be ~ 1 (always £ 1). Hence a sufficient 
criterion for fluorescence to provide negligible rotational excitation is for 


u e r J' , e f, J ,,B e ,, J n , e'J’ < <A J M , J"+l 
for all e', J'. 


(6- 3) 


B. Vibrational Fluores cence 

We will assume that the galactic radiation field near X ~ 4. 9/x 
(the wavelength of the CN v= 0 -* lvibrational transition) is predominantly 
dilute starlight. Speculation aboutThe existance of radiation from young gal- 
axies (Partridge and Peebles 1967) suggests an additional contribution which 
is at most the same order of magnitude, so we will neglect this component. 
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The energy density of starlight at this wavelength is given by 

—28 —1 -1 
u(4. 9 p.) = 8. 8 x 10 erg cm Hz 

The Einstein B coefficient for this transition is calculated in Appendix 3, and 
from it we calculate that 
R 


v=0, 1 


-12 -1 
< 4 x 10 sec 


This rate is totally negligible with respect to the excitation rate necessary to 

0 

produce the observed 2. 7°K excitation = 8. 7 x 10 . Thus vibrational 

excitation falls short by six orders of magnitude in accounting for the observed 
excitation. 


C. Electronic Fluorescence 


For electronic transitions in which one electron is involved, we can 

/ 

write the Einstein B coefficient in terms of the electronic oscillator strength 

or "f value" ( IX Herzberg 1959) . Thus the upper limit for the excitation 

by optical photons of frequency v= v . . is given by 

© v J , e v J 


R *c 

J", J"+l ~ 


u(v, n s 2 W 1^1* US? I* 


VT 




mhc v (2J"+1) 


(6. 4) 


where are the Honl - London factors, • and ^ are the 


J"J 

Franck - Condon factors. 

Then the usual sum rules apply 


££ I® 1 '™* 1 , | 2 = 2J" + X 

■3 f M'M ,r 


(6. 5) 
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D. 1^ 

V* 


(vib) 

v'v” 



1 , 


( 6 . 6 ) 


s w - 1 (6 - 7) 

e' 

Since u(v) is a slowly varying function of v, and the excitation will be pre- 
dominantly via the (0, 0) band, we may sum over allowed rotational and vibra- 
tional transitions to get 


R 


01 


u(v) 


2 , 

f ...» 


e e 1 


mhcv 


( 6 . 8 ) 


where 


v v e ' v=0 J"=l, e"v" =0 J" =0’ 


C a CN Band Systems 

CN has two band systems in the visible, the well known red 
2 2 2 2 

(X £ - A p) and violet (X 2 - B 2 ) systems. The f values for these 
transitions have been measured to be 0. 0034 (Jeunehomme, 1965) and 0. 027 
(Bennett and Dalby 1962) respectively. 

2 

Ultraviolet systems involving electronic states above B 2 have 

been studied by Douglas* and Routly (1954, 1955), and by Carroll (1956). Our 

\ 

current knowledge of the electronic term scheme of CN is shown in Figure 14. 

It is of considerable importance to us to know the oscillator strengths 
of the CN ultraviolet systems. Douglas (1968) has searched unsuccessfully for 
these transitions in absorption down to X — 1400A, in. a long path length 
"White-Cell". He has concluded that all of the oscillator strengths for these 
ultraviolet systems, which terminate m the ground electronic state, are less 
than the oscillator strength of the violet system (. 027). The CN dissociation 
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energy corresponds to \ ~ 1500 A (Carroll 1956), thus any state with a 
large oscillator strength for transitions from the ground state must lie beyond 
the dissociation limit. 

C, Upper Limits 

~b — ol 

It is revealing now to calculate the upper limit for R^. given by 
(6. 8) assuming f and a to be equal to 1. This is plotted in Figure 15 as a 
function of wavelength, for a radiation field corresponding to that of the galactic 
starlight background (Zimmerman 1964). For the red and violet bands the 
measured f values are also shown displacing these upper limits accordingly. 

Thus it is clear that for any band -system, the fluorescent excitation 

-9 -1 

rate due to the galactic background will be ~ 2 x 10 sec , or smaller. This 
is ~ 6000 times too small to yield the required excitation, and we may thus 
conclude, without further discussion, that .the galactic starlight background is 
not responsible for the CN excitation. This is a most important point, since 
this is the only process which we have been able to find which might account 
for the observed independence of T^CN) with respect to galactic latitude, 
longitude or spectral type of the star against which the CN is observed. 

■Now consider the situation in which the star is the source of fluores- 
cent radiation. Figure 15 also shows the upper limit for R given by (6. 8) 
for a radiation field corresponding to that of a line blanketed BOv star at a 
distance of lOpc (Hickok and Morton 1968). Again we assume f and a = 1, 
and see that at this distance, the starlight intensity is of the same order as the 
galactic starlight. Thus if the CN is not circumstellar, the radiation field 
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of the star used to observe it is of insufficient strength to account for the ob- 
served rotational excitation. 

In order to have the fluorescent excitation of the star play a dominant 

_5 

role in the rotational excitation of the CN, we need ~ 10 sec . There 

are two consequences of an excitation rate this high. First, the flux necessary 
to yield this high a rate would require the star to be as close as 0. 2pc to theCN. 
This proximity is totally contrary to our present ideas concerning the location 
of interstellar clouds. A review of our present understanding of their location 
will be deferred to Chapters 7 and 9, however, the main feature which emerges 
from this work is that the clouds are not circumstellar, but rather broad layers 
spanning, in several instances, more than one star. For example, the inter- 
stellar CN lines' of X Persei and' C Persei appear to originate in the same 
interstellar cloud. These clouds appear to have the atoms and the molecules 
well mixed, and only a moderate radiation field incident upon them (see for 
example, Herbig. 1968a). 

Secondly, since we know that all of the f values are small 
electronic states which lie below the dissociation limit, the sum rule (6. 7) 
implies the existance of at least one electronic state above the ionization limit 
with an f value approaching unity. If a molecule makes a transition to this 
state, it will have a certain probability for dissociation. It is difficult to estimate 
what this probability is, but we believe a value of ~ 1 % to be a conservative 
guess. 

Then, since’ the f value of the dissociating state is greatest, this 
state will produce the fluorescence, and the molecules will be destroyed at a 
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rate ~ 10 sec per CN molecule. (They will be excited at a rate 10 sec , 
and 1% of these will dissociate. ) In order to maintain an equilibrium population 
of CN, we must have a continuous production rate for these molecules of equal 
magnitude. 

The production of interstellar molecules is unfortunately not well 

understood (see for example Dieter and Goss 1966), and only with difficulty is it 

possible to get CH production rates comparable to CH dissociation rates (~ 10 ^ 
-1 

sec per molecule, for the CH existing in the interstellar radiation field). 
Assuming similar processes to be responsible for the formation of CN, it would 
appear, then, hard to understand even the existance of the CN. Hence it appears 
highly improbable that the observed excitation could be due to the fluorescence. 
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CHAPTER 7 

COLLISIONAL EXCITATION IN H I REGIONS 

j * 

In this Chapter we summarize briefly those properties of H I and 
H II regions which are pertinent to the problem at hand, and we review work 
devoted to the location and environment of the interstellar molecules. We then 

' ' * i 

proceed to consider in some detail the possible rotational excitation mechanisms 
of the CN if it resides in an HI region of the interstellar medium. 

A. Difference Between H I and H II Regions* 

i. - 

Almost 10% of the galactic mass is in gas clouds ~ 4pc in diameter 
which comprise about 10% of the volume of the galactic plane. These tend to 
be irregular m concentration and composition, but are composed mainly of atomic 
hydrogen, (although H may predominate at higher densities). This gas is divided 
into two distinct types of regions: H I and H II regions (Stromgren, 1939,1948). 

H I regions are neutral atomic hydrogen with a kinetic temperature ~125°K, and 

4 -3 -3 

hydrogen number density rarely exceeding 10 cm and usually ~ 1cm . 

H II regions, on the other hand, are the emission nebulae which 

surround hot, luminous O and B stars which have sufficient ultraviolet flux 

to fully ionize the surrounding atomic hydrogen. These regions, then, have a 
4 -3 

temperature ~10 °K and density ~10cm . Since O and B stars are com- 


* Eor a review of the interstellar medium see Dieter and Goss (1966) and Vol. 
VII of Stars and Stellar Systems: Nebulae and Interstellar Matter (Middlehurst 


and Aller, 1968) 
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paratively uncommon, and stars later than B 0 have only limited H II regions 
associated with them, ' H n regions only involve ~ 10% of the mass of the inter- 
stellar gas near the galactic plane. 

A well known but noteworthy feature of the interstellar medium is that 
the boundary between the two types of regions is relatively, sharp (Stromgren, 
1939, 1948). This is due to the fact that the radiant energy density beyond the 
Lyman limit falls off inversely as the square of- the distance from the star, where 
as the total opacity due to the non-ionized component of the gas increases exponen- 
tially with the path length to the star. These two effects', when included m the 
equations -‘for ionization equilibrium produce'a very sharp transition between H I 
and H II. 


B. Earlier Work on the Place of Origin of the Molecular Interstellar Lines 

Prior to Adams' work (1943 } 1949) Merrill (1942, 1946) had suggested 
that there might be a difference in the place of origin of the molecular absorption 
lines and those of Na I and Ca II,’ known to be interstellar in origin. The latter 
would presumably originate in' the "true interstellar clouds", far from the star 
against which the lines were observed, while the molecular lines might be produced 
in circumstellar clouds, and be affected by the stellar radiation field. 

To clarify this situation, Adams measured the velocity correlation 
between the H and K lines of Ca II and the molecular lines. He found that 
all of the evidence suggested a common origin for molecular and atomic lines 
in the same clouds, and in the same portions of the same clouds. He showed 


that the molecular lines, when present, almost invariably correspond to the 
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strongest components of -complex H and K lines Thus there is a close relation- 
ship between the Ca II and the molecular lines as regards both position and in 
most cases, relative intensities as well. 

Also contributing to his conclusion were the presence of the same 
molecular lines with nearly the-same intensities m wide double stars and numerous 
pairs of physically unrelated stars. He also found no evidence that the occurrence 
of the lines is related to the spectral types of stars which appear to he nearest 
the clouds of their origin. 

To further support this work, Bat.es and Spitzer (1951) subsequently 
showed that Adams'- visual estimates of- the intensities of the molecular lines 
yield a correlation with the intensities of the interstellar reddening, the latter 
is presumably due to the presence of grains. * Taking advantage of the new and 
more accurate atomic and molecular f values,- Herbig (1968a) has recently con- 
sidered this problem for one particular star f Ophiuchi. This was chosen because 
it has a single well resolved interstellar cloud m front of it He has arrived at 
a similar conclusion that the place of origin of both the molecular and atomic lines 
is an interstellar H I cloud. His arguements for this conclusion are as follows: 

* Bates and Spitzer (1951) also suggest that the stars m the Pleades represent 
a unique case, in that the CH+ m the spectra of these stars is probably circum- 
stellar. 
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1. The pattern of Ca IT absorption lines m £ Ophiuchi is repeated 
in. at least seven stars in this direction, many of them -well outside the H II 
region surrounding £ Ophiuchi. Since this is the only star in the group involved 
in an H II region, these lines must be formed in a rather extensive foreground 
layer not associated with these stars 

2. A similar arguement applies to H I 21cm emission feature 
which is not appreciably affected by, and extends well beyond the £ Ophiuchi H II 
region. This feature is identified with the same cloud containing the other absorb- 
ing molecules. 

3. If the Na I lines were formed in an H n region, the increased 
electron presure would depress the sodium ionization and imply an order of 
magnitude deficiency in the Na/H ratio. If one assumes an H I region as the 
origin of the Na I lines, one gets the solar system Na/H and K/H abundance 
ratios. 

4. A comparison of the 109 a and Ha emission lines originating 

in the H II region shows no turbulent width correlation with any of the interstellar 1 
species. 

Thus Herbig has concluded that for the case of £ Ophiuchi at least, 
which is probably the most important star for our work, the place of origin of 
the molecular lines is an interstellar H I region. 

C. Collisional Mechanisms in H I Regions 

We now turn to a consideration of the various collisional excitation 
mechanisms of the CN, assuming that it does reside in an HI region. Since 
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the kinetic temperature here is ~125?K, and particles -will have kinetic energies 

~4kT/rr = 229cm it is clear that we need to consider only direct rotational 

excitation.- Indirect rotational excitation by way of higher electronic and vibration* 

al states can be immediately ruled out -on energetic grounds alone, as, the thres- 

-1 -1 

holds for excitation of these states are 2024cm and '9114cm respectively. 

D. Excitation by H Atoms 

We -shall consider first the collisional excitation by neutral hydrogen 
atoms. We can write for the excitation rate 

E 01 “ -Kl (v >' 01 01 (v) dV ” N H v °01 “ N H r 01 (7 - J) 

where: 

N is the total H I number density 
H 

v is the mean atomic velocity 

<Tq^ is the total cross section for rotational excitation 

The calculation of m terms of molecular parameters is a difficult 
problem m the theory of atomic collisions. Fortunately, the problem is one which 
m another guise has been studied from both an experimental and theoretical point 
of view. This is the problem of microwave pressure broadening (Townes and 
Schawlow, 1955 p. 338). 

Virtually all of the experimental work of this subject has been devoted 
to an understanding of the interaction between chemically stable systems which 
can be conveniently studied under laboratory conditions, while H I and CN are 
both free radicals. The particular short range forces responsible for pressure 
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broadening, however, do not basically distinguish between free radicals and 

chemically stable molecules, and it is thus possible to estimate the cross section 

for CN rotational excitation by H to an accuracy of perhaps 50% by an appeal 

to the laboratory data of similar molecules. 

A rather similar process which has been described in the literature is 

that of the pressure broadening of the H Cl J = 0 -» 1 rotational line by collisions 

with He atoms at room temperature (Rank et al. 1963) . The resulting cross 

section for the J = 0 -» 1 rotational excitation has been calculated from the ex- 

-15 2 

penmental results to be 1. 6 x 10 cm (Gordon, 1966). 

H Cl, like CN, is a highly polar molecule with^ Tjri1 = 1 085D (Bell 

xi GL 

and Coop, 1938). The bond lengths of the. two molecules are also similar - 

r (HC1) = 1. 274A versus r (CN) = 1. 172 A (Herzberg, 1959, Table 39). 
e e 

Unfortunately the quantum-mechanical description of the collision 
between a polar molecule and a spherically symmetric atom results in a slowly 
converging expression for the excitation cross section (Anderson 1949). Since 
each term scales in a different way with the molecular dipole moment p , the 
atomic polarizability a, and the projectile velocity v, there is no precise way 
m which we can deduce the CN - H or CN - H excitation cross sections from 

A 

* 

those of H Cl - He. However on the basis of the leading term m the expression, 

due to the dipole - induced dipole interaction, we can estimate that an upper 

2 2/5 

limit, at least, for the cross section will scale as (an /v) 


^Actually for CN - H collisions there is a comperable term due to the direct 

interaction between the CN dipole moment and the H quadrupole moment. 

z 
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In addition to a Q for HC1 -■ He collisions, Gordon (1966) also gives 
the corresponding cross sections for OCS - He and OCS - Ar collisions The 
OCS electric dipole’ moment is small, and consequently we have not used this 
molecule as a comparison since interactions other than the dipole induced - dipole 
type probably are moreimportant. But it is encouraging to note that these data 
scale according to our prescription. 

The polarizability a of the colliding atom may be calculated from 
the quadratic stark effect of its ground state and has been calculated for both 
H and He, as well as having been measured for the latter atom (Bethe and Sal- 
peter, 1957). It is found that a^/a^ e = (9/4) / • 074 = 30. 4. 

The ratio of the velocities of 125°K H atoms to 300°K He atoms is 
1. 29. Thus we expect a cross section for this process of a (H - CN) = 

•ct (He-HCl) x 4/5 =7.1 x l(f 15 cm 2 

Taking Herbig's (1968a) upper limit of 900cm for the neutral hydro- 
gen number density in the cloud containing the CN in front of Q Ophiuchi, and 
_ ^ 

v = y8kT7(nm^)' = 1. 6 x 10 cm/sec for the mean velocity of 125°K hydrogen 

”6 — 1 

atoms, we get the excitation rate R « 1. 03 x 10 sec which is smaller 
~ -6 -1 

than B q ^ u^ q - 8. 7-xlO sec by a fact or of 8. 4. Thus it appears improbable , 

i 

although not by a large margin, that collisional excitation by hydrogen atoms could 
account for the observed rotational temperature . At most a small portion of the 
observed excitation can be attributed to collisions with H atoms * 

* There also exists the possibility that a fraction of the interstellar hydrogen may 

be in the form of H molecules The ratio of the polarizabilities of He and H 
£ 2 

is given by a / = 3. S3 (Margenau, 193S), so we expect a cross section 

a 0l (CN " H) / a oi (CN " H 2 } = L 97 ' 
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E . El ec tr on D ens ity 

In H I regions, the ultraviolet flux is of insufficient magnitude to 
ionize the hydrogen gas there, however, there is sufficient radiation to ionize 
such more loosely bound but considerably less abundant species as C, Mg, Si, 
and Ee, primarily C. It is generally assumed that the electron density is due to 
photoionization of these elements (Middlehurst and Aller, 1968). More recently, 
cosmic rays also have been suggested as a possible source of ionization (Potash, 
1968; and Spitzer and Tomas ko, 1968). 

-3 

In H I regions of typical density ( 1 cm ) the electron density is 

i 

thus quite small. For example recent dispersion measurements made with the 

-2 -3 

help of pulsars yield electron densities of order N g « 2 x 10 cm (Shuter et 

al, 1968) averaged over hundreds of thousands of parsecs in the galactic plane. 

However, in the cloud containing the CN in front of Q Ophiuchi where the density 

3 

may be as high as 900 atoms/cm , the density of these elements and hence the 
corresponding electron density will be proportionately higher. From a consider- 
ation of the expected carbon to hydrogen ratio, Herbig (1968a) has estimated that 

_3 

here it may be as high as N ~ . 4cm . 

6 

If either of the optical lines of Na II or Ca I were observed, one 
could calculate the actual electron density in this cloud from the equations for 


It thus appears that two H atoms are roughly four times as effective m exciting 
molecular rotation acting alone, as they are when bound together m a molecule. 
Thus we have assumed, for the sake of an upper limit, that all of the hydrogen 
is in the atomic form. 



75 


ionization equilibrium. Since only Na I and Ca II are seen (due to thei r differ- 
ence m ionization potentials), to use the Na i/Ca II ratio one would have to 
assume the abundance ratio of Na to Ca. However the calcium abundance is 
often anomalously low, - and in the case of the Q Ophiuchi cloud Herbig points 
out that this effect is probably exaggerated. The easiest way to measure the 
electron density here would be to attempt to observe the Ca I resonance line at 
4227 A for this star. Herbig has calculated from his upper limit of <3m? for 

_3 

the equivalent width of this line that < 1. 9cm 

It is unfortunate that this line was just off the bottom of the CH+ syn- 
thesis described in Chapter 5 of this thesis (this synthesis terminates at 4228A). 
Further effort should be made for the detection of this line. 

F Cross Sections for Excitation of CN by Electrons 

In order to estimate the excitation cross section for 125°& electrons, 
we must consider the problem quantum mechanically, and will use the Born approx- 
imation for this. ■ Because these -results are used later when we consider H II 
regions, we will solve this problem for the general case of a diatomic molecule 
with both spin and orbital angular momentum, and specialize to the case of a rigid 
rotor for CN. Although Massey- (1932) derived a similar expression earlier, 
his results appear to be m error by a factor of 2/3, hence, we shall present the 
calculation' here. 

The usual condition for the validity of tins approximation applied to 

' i ' t 

2 

the excitation of a point dipole is that 2fjiem/ft « 1 , or p « 1. 27D (Massey, 1932). 
While we adopt p(CN) = 1. 45D, and p(CH) = 1. 46D, (Table 3. 2) it is known. 
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however, that Massey’s condition is often too stringent a restriction on the use 
of the Born approximation. Altshuler (1957) for example, obtained good agree- 
ment with experiment by calculating the diffusion cross section of NH and H O 

o A 

in the adiabatic Born approximation, and these molecules have dipole moments 
(1. 468D and 1. 94D respectively) ,as large and larger than CN and CH. Mittleman 
and von Holdt (1957) have moreover calculated the exact cross section for elastic 
scattering from a point dipole, and find that even at [X = 1. 46D the diffusion cross 
s'ection is only 60% higher than the result in the Born approximation. 
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where 0 is the scattering angle of the electron , and we take the z axis to be m 
s 

the direction of 6k = k. - k as shown, in Figure 16b. Then the matrix element 

~ ~i t 

in (7. 4) may be expressed simply as 
i6k • R 


v a m : T / 

J_r 


ep j ^i6kRcos0 


< y - 1 etL£os_ x , A , ^2 


A K S JM 1 > 2 ! ^AK'SJ’M ; Rsm@d@d $ 


(7.5) 


sin © < ijr 


A KSJM 


° 0SX l*AESJ'M >,3Rd ® ,iS 


2 2 

with all other terms defined m the Figure. Note also that (6k)' =k + 

2 1 
k^ - 2k. k^ cos 0^'. -We can now expand cos \ m terms of spherical harmonics: 

. 1 

4tt v ' \ ix 

cosy-—- S (-D AV 8 * 1 (7.6) 

M - 1 


Since 5 appears nowhere else in (7. 5) the only term which will contribute to this 
sum is the one with p = 0. Hence we have for the matrix element, using the 
Wigner-Eckart theorm (Edmonds, 1968, equation 5.4.1), 


_ co -pr 

nr, _ 2ne p /X-Vo f C i6 kRcos© 

•K- “ q /4 tt/3 J J e sm© cos©d 0dR 


if L" 0 0 

^ /V 

'<AKSJM| Y l AK* SJ'M' ) 


(7. 7) 


= (-D 


J-M 4nep i 
L 3 6R 


J 1 J' 
-M 0 M' 


< AKSJfl C(l) II AK'SJ' > 


Using Edmonds (3. 7. 8) and integrating,equation (7. 4) then reduces to 


c-2 2 2 a a „ k ti sin0 dB 

a KJ,K'J’’ V I <AKSJ|!0(1) 11 AK>SJ'>I j4j % S 

3(2J+1) A J l 0 k + k -2k k cose 

if if s 


„ 2 2 2 
OTTin e p 


3(2J+1) h 


1 AKSJ |l C(l)|l AK'SJ) j . - In 


V fc e 

k-k ( 


(7. 8) 
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The evaluation of the reduced molecular matrix element is carried 
out in Appendix 2, hence for a diatomic molecule, such as CH, intermediate 
between Hund's cases (a) and (b) 


r A ^ = vcr A * = MIL A 

KJjK’J 1 U KJ,K'J' n 2 

3h v 



(2J'+ 1) 


• | K 4k (2K'«) 1/2 (2K + 1) 1/2 


K' J 1 s|(k 1 K'\ 

J K ljl-A 0- A y 



(7. 9) 


where k ~mv/ft, and k. = k. - 2mtt> '' " /A. When specialized to the linear 

1 II KJ ,iv3 

rotor (7. 9) yields 


JJ 


3h% 


2 2 k+ k 

i f 

k.-k 

l f 


. Sum e , i 

. " o 1 


<2J'+1) 


J 1 J*' 
0 0 0 


(7. 10) 


which is a factor 2/3 smaller than Massey's (1932) result. 

Using the approximation that 

(k. + k ) / (k. - k ) r s 2mv 2 / Auj. 
i f i f i] 

we can average the cross sections over the Maxwell-Boltzmann velocity dis- 


(7. 11) 


tnbution. Thus for a molecule such as CH 

r 


(r A ) - - 

v KJ, K'J 1 ' 


32 y 2 e 2 


3hv (2J+1) 


In 


-2 
tt mv 


2hai 


I" Y 


K'J',KJ 


S U K'K' 4k < W 7 W)* 7 * jK* J' sjl' K 1 K ,X 
and for the rigid rotor 


t2 


(7. 12) 


J K 1 -A 0 A 


2 2 

(r > = 

JJt 3 A 2 v 


In 


tr mv 


2 *“k'j’,kj 


(7 13) 


' where y-0. 577 (Euler's constant) 
and v = ,/ 8kT /ran 

The use of the approximation (7. 11) results in an error of less than one percent 
in the average rate, and we neglect the slight difference between the electron 
' mass and the reduced electron mass. 
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G. Excitation Rate from Electrons 

Application of equation (7. 13) for Herbig's (1968a) estimate of 

-3 -7 _i* 

N 4cm yields the excitation rate R.. »i5. 6 x 10 sec , which is 15 
e 01 

times lower than B^u , On the other hand it is only 3. 3 times smaller if 

-3 

we use Herbig's upper limit of N < 1. 9 cm . Thus it appears improbable 

G 

although again not by a large margin that slow electrons could provide sufficient 
excitation to account for the observed rotational excitation of CN . It also appears 
that for a normal C/H abundance ratio, electrons are approximately as effective 
as H atoms in exciting molecular rotation. 

Since the Born approximation satisfies electric-dipole selection 
rules, the rate equations may be solved by an equivalent increase of radiative 
intensity at the corresponding wavelength (this proceedure was described in 
Chapter 2). 

The increase in temperature due to excitation by slow electrons 
is given in Figure 17 as a function of electron density and observed temperature, 

I 

and for H I atoms assuming that the CN - H interaction predominantly obeys 

_3 

electric-dipole selection rules. For Herbig's estimate of N = . 4cm , the 

© 

correction to T and T due to electron excitation is 0. 14°K and 0. 12"®K 

XU 

respectively, just comparable with the statistical error of the temperature 
measurement. 

_3 

If the electron density is as high as Herbig's upper limit N <1. 9cm , 

e 

the corrections are still less than 0. 75°K and 0. 77° K, with the result that the 
major part of the excitation must be due to some other process. Further study 
of x 4227 A of Ca I should allow the upper limit to be considerably depressed. 
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H. Comparison of with. to Determine Electron Density 

Since collisional excitation is a nonthermal process, the resulting 
distribution, of level populations will not be thermal. In principle, this departure 
from thermal distribution can be measured and additional information about the 
collisional mechanism obtained. 

The expected T for excitation by a mixture of microwave back- 
ground and 125° K electrons has been plotted in Figure 18 as a function of 
electron or H I density and observed T It can be seen that the departure 
from thermal equilibrium of the J = 2 level population is too small to be measured 
with the present precision of T measurement; the temperature differences 

ui 

are at most « 0. 1°K for reasonable electron densities, comparable to the 
statistical errors involved. Thus with the present observational accuracy, a 
partial contribution from electrons or H atoms cannot be ruled out 


I. Excitation by Ions 

In the next Chapter, we will describe a numerical calculation of the 

cross section for rotational excitation by particles with the mass of a proton or 
\ 

greater. Figure 19 shows the results of this method applied to the cross section 

for excitation of the CN J = 0 -» 1 transition. In this curve the abcissa is in 

4 

units of the most probable velocity of 10 °K protons. For an H I region the 
excitation will be primarily due to C II ions at 125°K, hence the most probable 
velocity of these will be down by a factor 1/31. 
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It can be seen that the cross section is falling off very rapidly with velocity, 

and at the most probable velocity of the ions, it will be at least two orders 

4 

of magnitude below the value at the most probable velocity of 10 °K protons. 

In the next Chapter it will be shown that a number density of approx- 
3 4 

imately 1. 37 protons/cm at 10 °K is sufficient to provide 3°K excitation of 

the CN. Assuming charge neutrality in the H I region the ion density will be 

the same as the electron density ( 4e/cm 3 ). Thus the excitation rate here 

4 

will be smaller by a factor of approximately 10, and will clearly be negligible. 

J. Discussion 

It is thus seen that if the CN resides in an HI region typical of the 
major interstellar cloud in front of Q Ophiuchi, its excitation probably provides 
a good measure of the background radiation at 2. 64mm. No other excitation 
mechanisms appear to be of sufficient strength to provide this excitation, al- 
though the margin for e and H excitation in the case of Q Ophiuchi is not as 
high as one might like. There is good evidence that the CN does m fact reside 
in an HI region, as we discussed earlier in this Chapter. Further discussion 
of this point will follow in Chapter 9 after we consider the alternative possibility 
of the CN residing in an HH region. 
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CHAPTER 8 

COLLISTONAL EXCITATION IN H II REGIONS 

111 H II regions virtually all of the collisions experienced by a 
molecule will be with charged particles. We will now show that, as a result 
of the long range of the Coulomb force, the excitation cross section is sufficient- 
ly large that all of the observed CN excitation could be produced in an H II 

_3 

region of reasonable ( ~ 1cm ) density. 

Our analysis of Chapter 2 showed that even though the CN may not 
yield reliable temperature measurements, in this case it will set a reliable 
upper limit to it. In addition we will now show that if protons are the sole ex- 
citing agent for the interstellar CN, the resulting rotational temperature will 
be a sensitive function of the density of exciting protons. Since the interstellar 
medium is a highly inhomogeneous place with respect to density, the observed 
invariance of the rotational temperature against eleven stars (see Chapter 9) 
speaks strongly against the suggestion that the excitation is collisional. 

We will calculate here the CN J = 1 and, J = 2 and the CH KJ = 

13/2 level populations which result if the CN J = 1 population is due solely 
4 

to excitation by 10 K protons. Following Chapter 7, we are interested to see 
if the present observational upper-limits on these levels might be useful in 
shedding further light on the CN excitation mechanism. Hence our calculations 
will be for the most complicated case of interest, a diatomic molecule, such 
as CH, with both spin and orbital angular momentum which lies intermediate 
between the Hund's case (a) and (b) coupling schemes. We can easily specialize 
this to the simpler coupling scheme of CN. 
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A. Semi Classical Approximation 

Our calculation of the excitation of molecular rotation will treat 

the projectile trajectory classically. Such a proceedure is usually valid in 

scattering theory when the impact parameter b is large with respect to the 

' projectile de Broglie wavelength )*= h / Mv (here M and v are the projectile 

0 

mass a velocity relative to the scattering center). For v « 10 cm/sec, the 

maximum excitation of CN and CH occurs when b «50A, so this condition is 

well satisfied. In addition, since the projectile energy is ~lev, while the lower 

-2 -3 

energy level separations of CN and CH are ~ 10 to 10 ev, the energy 
loss of the projectile will be small. Hence its deflection will also be small, 
and we may assume it to follow a classical straight line trajectory. 

The transition rate R for a density of N incident projectile 

Jd P 

(1ST (v) at velocity v>) is given by 

P , 

00 

R JJ' = Vjj’ = *( VCT JJ. W N p ( v ) (8. 1) 

where v is the relative velocity of the molecule and the projectiles, and 

^TJ'^ is the (velocity dependent) cross section for excitation of a molecule 
mitally in level J to level J’. The cross section is found by summing the 
transition probability over all initial and final states, integrating 

over impact parameter, and dividing by the statistical weight of the initial 
level. Written in the differential form this is 

23 + 1 d y 

2,0 db M’M “ 

The dynamical variables of the molecule, on the other hand, must 




be treated quantum mechanically. The molecular state vector during the 
encounter with the projectile is described by Schrodinger's equation m the 
form (see for example Schiff 1955, p. 196). 
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l h 


da 

l 

dt 


£ 


iyj. t 

a j 3C (t) e 13 


(8. 3) 


Even with the simplifying assumptions already made, no analytic 
solution exists for these coupled differential equations for arbitrary values of 
v, b and the molecular constants B or |t We must therefore solve these equations 
numerically, and use a Runge Kutta scheme to do this (see below). First the re- 
quisite matrix elements of the interaction Hamiltonian JC'(t) must be calculated. 
The excitation cross sections may then be found numerically, and averaged over 

A 

the velocity distribution of the 10 protons typical of an H II region. Since 
the computer time required for the Runge Kutta integration is large, for a check 
of the resulting rates at other temperatures we use a hybrid of first-order per- 
turbation theory and the quantum mechanical sudden approximation. This has 

4 

yielded the same cross sections to an accuracy of ~ 1 % near 10 °K (Thaddeus 
and Clauser, to be published). 


B. Matrix Elements of Interaction Hamiltonian 

The matrix elements of the Hamiltonian are taken to be those of the 
interaction of a point charge with a point dipole. Following equations (7.3) and 
(7. 6) we take the operator to be given by 


3(7 (t) ~ 3R^ ^ Y (®> s ) 

H=-l ^ w 


(8,4) 
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where JR, ® , and § describe the time dependent spherical polar coordinates of 
the position of the projectile in the laboratory frame with origin at the molecular 
center of mass, and 0 and (p represent the spherical polar coordinates specifying 
the molecular orientation in the laboratory frame. 

For a classical straight line trajectory of a projectile with impact 
parameter b and relative velocity v, the Y (®(t), $(t), ) may be written for the 
coordinate scheme shown in Figure 20. as 


1/2 


1 ± 


, (t) = + (3/8 t t) ■ vt 

1 (v^ + b 2 ) 


1/2 


Y 10 ( t} =(3/4? r) 


1/2 


(vV + b 2 ) 


1/2 ( 8 * 5 ) 


Hence the matrix elements of the Hamiltonian (8. 4) may be written 


^ * l/2 

< AKSJM i 31' (t) I AK' SJ' M' > = (4?r/3) jie. 

R ? (t) 


2) (-l)^Y (t) (AKSJMI C (0, 0)1 AK' SJ' M' ) 

_ Li 


( 8 .(?) 


p=-l 


The ^ are the components of a tensor operator, so we can apply 
the Wigner-Eckart theorem (Edmonds 1968, equation 5. 4. 1.) to factor off the 
M 1 , M dependence. Thus 

<A KSJMI y' (t)i AK' SJ' M' > = ( 4^/3 )^ 2 _ye 

R 2 (t) 


(8.7) 


<AKSJflC(l)l!AK'SJl> S Y, (t) (~1) J M + ^ ( J 1 J ' 

« pi=-i l!i v M 


The reduced matrix element in (8. 7) has been evaluated in Appendix 


2 so that the matrix elements of 30 may be written 
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- - . 1/2 , . A +S+ J +J -M+1 /ot , .-.1/2 9T ..1/2 

(AKSJM I 3C'(t) ! AK'SJ' M') = (4 tt/3) (-1) ‘ ( 2J +1) ( 2J !) 


$» Seif V t; l-M-UM'j 

. v u =- 1 


( 8 . 8 ) 


T U- n*. , 2K . ..l/a 1,1/2 K J' s / K 1 K' 

K ^ K U K’K U KK (2K ’ J (2K+1) j J K lj \-A 0 A 


Letting A = S = 0, K* = J 1 , K = J, and using Edmonds (1968) 
equation 6. 3, 2), we obtain the matrix elements for the linear rotor: 

<JM i y' (t) | J' M ) (47r/3) l/2 (-1) M (2J' + l) l/2 (2J + l) l/2 


J 1 J 
0 0 0 


Me 


£ (-if y 


R^(t) M = "l 


J 1 J' 
lM V-M -M M' 


(8.9) 


C. First Order Perturbation Theory 

The differential cross section of equation (8. 2) for a molecule such 
as CH has been calculated using first order perturbation theory* to be (Thaddues 


and Claus er, to be published) 

dr KJ, K'J' v dQ KJ,K' J' = e 0>fr j, 

db ».2 3 


db 


(2J + 1) 


3ft v 


jp K u k'K' W !k (2K ' + 1) 1/2 (2K + 1) 

-0 


1/2 jK' S' S | /K 1 K'' 

/J K 1 Uoa 


n 2 


( 8 . 10 ) 


*It should be noted that first order perturbation theory is in all respects equivalent 
to analyzing the pulse associated with the passage of the projectile into virtual quanta. 
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where K Q and" are modified Bessel functions. The transition probability 


for this process obeys the usual symmetric-top selection rules for electric 

* 

dipole transitions. To apply equation (8. 10) to the' simple rotor we again let 
A = S = 0, K 1 = J 1 , and K = J, and use Edmonds (1968 equation 3. 3. 2) to eval- 
uate the 6j symbol, 

dr jj’ da jj' 8 ? 2 e ^ , 2T+1 v/j 1 J't 
db V db 2 3 ( j x) \0 0 0 J 

3fi v V . / (8. 11) 




2 rj-z 


Again, the 3j symbol leads to the selection rule AJ = ±1 for the linear rotor. 

Figure 21 shows the application of equation (8. 11) to the J = 0 -* 1 
transition of CN, and Figure 22 shows the application of equation *(8. 12; to the 
KJ = 1 1/2 1 3/2 transition of CH. In both cases a logarithmic divergence 

of the cross section results at small impact parameters. 


D. Quantum Mechanical Sudden Approximation 

The integrated cross section has been calculated for the collision 
of a proton with a molecule such as CH using the quantum mechanical sudden 
approximation, and the result is 


r KJ, K'J 


° KJ, K J' 


2 2 
TT ft V 

3tie 


(2J'+1) 


£ (2L+1) / 
L 0 


2 

J L+l/2 * bV 


b 2 db 


£ U" u! (2K'+1) 1/2 (2K+1) 1 ' /2 
K’K K'K’ KK t 


| K 1 J 1 S | / K L K 1 
l J K L j l- A .0 A 


( 8 . 12 ) 
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The integral over impact parameter in equation (8. 12) (Grobuer and Hofreiter 
1958, p. 202) converges for L > 1. If K and K' are such that the second line 
of {8. 12) vanishes for L ^1, we can write 


2 3 


KJ 


> K'J' - <»- «, E ,-f^i 


„ 2 
3fi v 


L 


(L-l)L (L+l)(L+2) 


(8. 13) 


£ U K-k U KK ( 2KT+1 ) 1/2 (2K+1) ]k- J' S l ( K L 


K'K 


J K L V-A 


L K'\ 
0 A y 


For this case, letting A = S = 0, K' = J', K = J we obtain the excitation rate 
for the simple rotor 


_ _ 8i yV T 12L + il_ h L J'\ 2 

JJ' _ 2 X> f' (L-l)L(L+l)(L+2) \0 0 0 ) 

on v .L 

On the other hand, if the second line of (8. 12) does not vanish for 
L ^ 1, the rate diverges at large b, and we cannot evaluate the integral in 


(8. 14) 


(8. 12) from b = 0 to ■» we use the latter approximation. 

Using (8. 14) the sum of the many-quanta transition rates from a 
given level to all other levels can be calculated exactly. For excitation of a 
rigid rotor from the ground state to a level with J s 2 (8. 14) yields 


8TTM 2 e 2 (2J + 1) 

3fi 2 v (J - 1) J (J+l)(J+2) 


(8. 15) 


and by means of partial fractions we may write 


CO 


£ 

J=2 


OJ 


_ 8hA 2 

3* 2 v 


CO ' 


3/2 £ 
J,=2 




(8. 16) 
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The series terminates and is .easily evaluated:’- 


„ 2 2 

8ttM e 4 _1_ 
J=2 ° J 3ft 2 v 3 


E 


In the same way it is found that 

CO 

E 


and 


(8. 17) 


_ 8 n/£jjL_ , JL 

1 -I T 2 

J=3 3fi v 40 


{8. 18) 


^ 2 2 
E r 2J = 8m e _J_ 

J=4 3 hv 21 


(8. 19) 


Finally, the simple l/v velocity dependence of the excitation rate 

in the sudden approximation allows the Maxwell-Boltzmann velocity average to 

be done analytically. The average rates { r_l * ) and (r ) are given 

Xvd J J J 

2 2 2 

by equations (8. 13) and (8. 14) with the factor '8tt fJ- e / 3fi v replaced by 


0 22 
32p. e 

n 2-< 

3 h v 


where the average reduced projectile velocity 


( 8 . 20 ) 


l/2 

v = (SkT/ir m*) (8.21) 

and m* is the reduced projectile mass, and T is- the kinetic temperature of the 
projectiles (t emperature of the H II region). 

Figures 21 and 22 show the differential cross section in the sudden 
approximation, the derivative of equation (8. 12) evaluated for the specific case 
of excitation from the CN and CH ground states. The divergence of the sudden 
approximation for AJ = ± 1 transitions (L £ 1) at large impact parameters is 
evident, but the finite value of the many -quanta cross sections, as revealed m 
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equations (8. 13) and (8. 14), is also illustrated for the J = 0 -* 2 CN transitions. 


E. Num e rical Solution of Schrodinger's Equation 

The descrete points shown in Figures 21 and 22 result from a num- 
erical integration of Schrodinger's equation (8. 3), and at least to the accuracy in- 
herent m the numerical analysis, they represent an exact solution of the exci- 
tation problem. 

The amplitude of the state vector in equation (8. 3) is first separated 
into real and imaginary parts, 

a (w) = p (w) + iq.(w) 


and a change in independent variables is made (Winther and de Boer 1965): 

sinh w = vt/b (8. 23) 

Since the matrix elements of the Coulomb interaction, given by equations (8. 8) 


or (8. 9) are real, Schrodinger's equation can then be written as an infinite system 


of first order coupled differential equations m the real variable w, 

id .b sinh w 


dq,(w) , 

J , cosh w , .i .. . . i.. 

= : Zu <3 (w) 0 


dw 


ftv 


q^(w) sin— 


dp (w) . 

b =25lUL (W)U> 


dw 


hv 


u> b sinh w 


q.(w) cos ~^- 
1 v 


to b sinh w 

* - P (w) cos-^ 

1 v 

(8. 24) 

«) .b smh w 

+ p.(w) sin-^ 

i v 


Here j and i denote the set of quantum numbers AK'SJ'M and A KSJM, respec- 


tively (orJ'M' and JM for the simple rotor), and the initial conditions are pre- 
scribed at a.( -m ) while the transition probability is calculated from a (+«>). 

] ^ 

In practice, the infinite set of equations (8. 24) are truncated at some 
rotational level of the molecule, such that the inclusion of higher levels has a 
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negligible result on the transition probability of interest. In addition, the inte- 
gration is confined to a finite range, -w Q £w rSw^. The resulting finite set 
of 2N real first order differential equations is then solved over a finite inte- 
gration range by the Runge-Kutta method, subject only to the initial condition 

q i(~ W ()) = 0 3X1(1 = 6 ik 5 or q i (" w 0 ) = 6 ik and P i< _v V = °’ f ° r al1 1 

(i. e. , the molecule is m a definite state k prior to the collision). The tran- 
sition probability due to the collision is then 

(w Q ) + q 2 (w Q ) (8. 25) 

The specific numerical prediction scheme was a generalization of 
Hildebrand's (1956, p. 236) equations(6. 16 7) through (6. 16 10) to the solution 
of 2N simultaneous differential equationa Writing equations (8. 24) in the form 

dp dq 

-jJ- - F , «vv w) aBd dw = G J< ! W W) (8 26 

and defining 



P. 

jn 



P.(w ) + iq.(w ) 
1 n J n 


(8. 27) 


then the equations for the a. , in terms of the a are the following: 

■J.n+1 ' jn 

P . =P + 1/6 (k +2 k +2 k +k . \ 
j,n+l jn ljn 2jn 3jn 4jn ' 


and 


(8. 23) 




j,n+l 


= q + 1/6 (m +2 m +2nu + rn ) 
V 1 ljn 2jn 3jn 4jrf 



where 
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and 


k n = 6 w F (p +k . /2, q +nti /2, w +h/2) 
3jn j in 2in H in 2in n 

k = 61 F (p. +k„. , q, +m„. , w +h) 

4jn j ' m 3 in in 3in n ' 


m, = 6 w G.(p. , q ,w ) 

Ijn j m in n 

rri = 5wG(p +k /2, q. +m_. /2, w + h/2) 

2] n ] m lm m Tm n 

m = 6 w G (p +k. /2, q +in /2, w + h/2) 

3] n j m 2m in 2in n 

rn . = 8w G (p +k , q +m„ , w +h) 

4jn ] 'in 3m in 3m n ' 


and 


6 w = the increment in w 

To calculate the cross section of excitation rate for a given initial 
level via equations ( 8 . 1) and ( 8 . 2), the Runge-Kutta integration must be repeated 
starting from each state of the level in turn, unless it is possible to exploit 
the symmetry of the interaction For the simple rotor, for example, the rates 
are symmectrical witn respect to ± M. 

For CH, where only excitation from the ground state is important, 
the Runge-Kutta calculation included the interaction between all the states of the 
lowest seven rotation-fine structure levels of the molecule, while for CN, where 
excitation from the three lowesi levels is of interest, the lowest ten rotational 
levels were included in the calculation. The integration over w was confined 
to the range -4 ^ w < 4, and the increment 5w m the Runge-Kutta scheme 
was adjusted in the range 0. 005 to 0. 1, depending on b, to provide the desired 
convergence of the solution. Adequate convergence was furthur demonstrated 
by varying the number of levels, and w , and by comparing the solution with 
the sudden approximation and first order perturbation theory m the appropriate 
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limiting cases. Using this scheme each of the points shown in Figures 21 and 
22 required approximately one minute of computation on an IBM 360/75 computer. 

Figure 21 shows for CN how the solution calculated from first order 
perturbation theory converges for large impact parameters to the numerical 
solution and for small impact parameters how the sudden approximation converges 
to the numerical solution. From a consideration of the orders of magnitude 
involved (Thaddeus and Clauser, to'be published) it might be expected that the 
hybrid technique would yield accurate results when applied to CN,but be of only 
marginal value for CH. This contention is in fact born out as can be seen by 
comparing Figures 21 and 22. 

Figure 21 also shows, specifically for the J = 0 -* 2 transition of 
CN how well the sudden approximation alone is able to describe many-quanta 
transitions. Comparable close agreement between the exact and sudden approx- 
imation solutions is obtained for the J = 0 — *3 and higher many-quanta CN tran- 
sitions not shown in the Figure. 
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F. Numerical Results 

We now turn to the excitation of interstellar CN and CH, and 
specifically calculate (a) the number density of protons and electrons necessary 
to provide the observed 3°K excitation temperature of the J = 1 CN level, and 
the temperature of the J = 2 and J = 3 levels which then results, and (b) the 
excitation temperature of the KJ = 1 3/2 state of CH which this particle density 
produces. Charge neutrality and a common kinetic temperature for protons and - 
electrons are assumed, and it is supposed that CN and CH are well mixed to- 
gether, and thus subject to the same density of charged particles. 

a; CN 


Since for this molecule the population of several levels is appreciable, 
we are obliged to solve the rate equations (2. 14) which reduce to 
dn T » 


= 0 = -n T (N J) R tti +A t _ ,)+ n A + N Yj n Tt R TlT 
dt J P jJZq JJ J, J-l J+l J+l, J p jT q J' J'J 


(8. 29) 


where n^ is the population of level J, N^ is the proton for electron) number 

density, and R is now the sum of the electron and proton collision rates. 

J J 

These equations are to be solved for N , n and n under the condition n /n 

P ^ O -L U 

= 3 exp (-2hB/kT ), B = □. 7 x lO^sec \ T = 3°K, and the normalization 
B B 

condition T} n =1. That is, we assume that when levels with J > 3 are 
J=0 J 


collisionally excited they immediately relax radiatively to the J = 3 level, and 
that the population of these higher levels is effectively zero. In. view of the rapid 
increase of the spontaneous emission probability with J this approximation pro- 



duces a negligible error m the final results. The rate for many-quanta 
transitions from the levels J = 0, 1, 2 to higher levels is thus simply added 
on to the transition rate to the J = 3 level. The last column of the rate matrix 
m equation (.8. 29) is thus 

* “ OD 

B ,r 4 j 3 >+ t £ (8 - 30) 

for J = 0, 1, 2, where (e) denotes electrons, and (p) protons. 

The electron rates were already calculated in Chapter 7 using the 
Born approximation. The summation does not extend over the electron rate 
since only A J = ± 1 transitions are permitted in the first Born approximation. 

Note that (r^?, ) and {r^> satisfy the principle of detailed balance, but for 
terms which involve the J' = 3 column, R does not. 

Jti 

The proton contribution to the rate matrix is found by two techniques. 

In the first, or approximate, method the many-quanta rates are calculated from 
the sudden approximation, employing equations (8. 14), (8. 17), (8. 18) and 

(8. 19), since the integral over velocity as given by (8. 20) is expressed in closed 

4 

form. These terms are listed m Table 5 for an H II temperature of 10 °K. 

For the AJ = ±1 (L ^1) approximate rates we employ the hybrid 

solution described in the section on the sudden approximation. Fine intervals 

for both the integrations over b (6b = 1A) and v(6v = 0. lv = 0. l(2kT/m*) 

m p 

were used where v^ is the most probable proton velocity), and no appreciable 
error is introduced by the numerical integration. The approximate single-quanta 
rates for CN , were calculated m this way for H II temperatures from 5900° 
to 15, 000°K, at 1000° intervals; the rates for 10, 000°are listed in Table 5. 
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Iii the second, or exact, method both the -single and many-quanta 
transition rates were calculated by means of the numerical integration of 
Schrodmger's equation. Since the calculation of dr/db at a given value of b 
and v requires a substantial amount of computer time by this technique (at 
least one minute), an economical scheme for the double integration must be 
employed. Simpson's rule was used for the integration over b, with interval 
6b ranging from 2A to 10ft. The second integration over the velocity distri- 
bution was then done by a four point quadrature scheme — a modification of 
Gauss-Hermite quadrature to deal with integration over the domain 0 to 
Since m terms of v^ the average over the velocity distrbution is 


= 4 tt 


■1/2 -3 r 2.2 2 

J exp (-v /v m ) r j Jf (v)v dv 


m 


0 


(8. 31) 


an n point quadrature scheme will exactly integrate the first 2n - 1 = 7 terms 

2 

in a power series expansion of r ' (v) v m terms of v. In the sudden approx- 

J J 

2 

imation r (v) v is exactly proportional to v, thus high accuracy might be 
dd 

expected with a small number of quadrature divisions. The divisions of v/v^ and 
corresponding weights used here are those calculated by Huang and Giddens 
>(1968). The adequacy of only four points in this integration can be demonstrated 
by comparing the quadrature velocity average of the approximate rates with the 
result of the fine-interval numerical integration given in Table 5: agreement to 
better than one per cent is found m all instances. 

The exact rates averaged in this manner over the velocity distribution, 


4 

for an II n temperature of 10 °K, are also listed in Table 5. The tendency of 
the average exact rates to be smaller than the average approximate rates is due 
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to the divergence of the sudden approximation at low velocities, where the exact 

rates instead vanish (since the low frequency photons associated with the passage 

of the projectile cannot produce transitions). 

Using'the results of Table 5 and the molecular constants* of Table 

2, the rate equations (8. 29) truncated at J = 3, are solved numerically, and 

yield Figures 23, 24, and 25. First of all, we find from the exact rates that a 

-3 -3 

number density = 1. 37cm (and 1. 37cm ) at a kinetic temperature of 
4 

10 °K is adequate to account for the observed 3°K excitation of the J = 1 CN 
level. Such a density of charged particles is easily attainable in an H II region, 
so that the observed 'excitation could be produced collisionallyaf the CN resides 
in such a region. 

Next, in Figure 24, we see that if collisions excite the J = 1 CN 
level to 3°K, the J = 2 level will show very nearly the same temperature, again 
over a wide range of H II temperatures. The first significant departure from 
thermal equilibrium occurs with the J = 3 level, whose fractional population 

_3 

is only ~ 10 , and as such is not observable with present techniques. 

Finally, Figure 25 shows that the temperature T of the J = 1 CN 
level is a rapidly varying function of N . Since the density of the interstellar 
medium is highly inhomogeneous, the invariance of T Q against eleven widely 
separated stars would tend to disprove the hypothesis that the excitation is 

* Note that the calculations in this Chapter used the value jJL~l. ID (Arpigney 1964) 
for the CN dipole moment instead of the more recently measured value 1. 45D 
(Thomson and Dalby 1958). The reason is that this was the best available value 
at the time the work fur this Chapter was done, but the slight difference between 
these values should not affect appreciably the final values for the proton density 


and T which we calculate. 

jUJ- 
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colli sionally induced, by heavy charged particles. 

b) -CH 

For this molecule we are only interested in the temperature of the 
KJ = 1 3/2' first excited state. The numerical integration of Schrodinger's 
equation, averaged over the velocity distribution by quadrature, 

4 

yields for the proton excitation rate at an H II temperature of 10 °K, 

r ^/2 ^ 3 ^2 = = 9. 97 x 10 8 cm 3 /sec. The many quanta transitions to 

higher levels, and the single quantum transition KJ =■ 1 1/2 - 2 3/2, produce 

only slight additional excitation of the KJ = 1 3/2 level, and will be neglected. 

For the electron excitation rate we use equation (7. 12) and find 

r^^ ^ ~ = 4. 11 x 10 8 cm 3 /sec. The total excitation rate. is thus 

R = r^ + r^ = i. 41 x 10 ^cmVsec. Since the Einstein coefficient for the 

-5 -1 

KJ = 1 3/2 level is 9. 184 x 10 sec , the fractional population of this level is 

n = N R/A = 2. 10 x 10 _a (8. 32) 

P 

_g 

for N = 1. 37 cm , and the excitation temperature in turn is 
P 

T = hv/kln (2/n) = 3. 76°K (8. 33) 

11 

for v = 5. 336 x 10 Hz. The resultant fractional population is, however, still 
well below the limit of the present threshold of detectability. 
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• CHAPTER 9 . 

FURTHER EVIDENCE THAT THE CN EXCITATION ORIGINATES 
NON-COLLISIONALLY IN AN H I REGION 

In the previous three Chapters we have seen that neither infrared 
nor optical photons nor collisions with atoms, ions, nor electrons are sufficient 
to account for the observed rotational excitation of the CN, man HI region. 

We found, how ever, that if the CNis man HD region, the observed excitation 
could be produced solely by collisions with charged particles. The location of 
the molecules is thus of crucial importance to our arguement. 

As we learned in Chapter 7,- the question of the place of origin of 
the molecular lines was first investigated by Merrill (1942, 1948), and subse- 
quently by Adams (1949) and Bates and Spitzer (1951). Herbig (1968a), has more 
recently considered this question, with specific reference to the interstellar lines 
of Q Ophiuchi.With the possible exception of the CH and CH+ interstellar lines 
m the direction of the Pleades, all of the evidence pointed to an HI origin of tne 
molecular lines. 

We will now present further evidence that the CN excitation originates 
non-collisionally in an HI region: 

1 Probably the strongest evidence yet that the excitation is not 
collisional is the observed invariance of the CN rotational temperature in eleven 
stars. 

2. We present the results of H I 21cm observations made m the 
direction of these stars. The CN velocity profiles are shown to correlate well 
with the H I velocity profiles in the same direction. This is consistent with the 
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contention that the CN is well embedded in an , H I region of the interstellar 
medium. 

A. T ^ Invariance 

Let us first consider the observed invariance of the rotational temp- 
erature for eleven widely separated stars. This is probably the most cogent 
arguement that the place of origin of the molecular lines is not an H II region, 
and that the source of the rotational excitation of the CN is not collisional. 

Listed in Table 4 are the optical depths of the interstellar CN lines, 
as well as the rotational temperatures calculated from them. The saturation 
corrections used in calculating these values were described in Chapter 3, and 
appear m this Table under the column labeled AT. 

In Chapter 8, it was shown that the assumption that the excitation is 
4 

due solely to 10 °K protons (the dominant exciting agent in an H II region) yields 
a rotational temperature which is for example, from 2. 5° K to 3. 5° K as the 

_3 

proton number density varies from only 0. 9 to 1. 9cm (see Figure 24). 

Considering that the density of the interstellar medium is highly 
inhomogeneous, one would expect a wide variation in CN rotational temperature 
from one cloud to the next if its excitation were local. However, the present 
data do not support this suggestion. On the contrary, the observed CN rotational 
temperature appears to be invariant and independent of cloud location. For all 
of the entries m Table 4, including a total of eleven stars, the data may be 
said to be consistent with the idea that the spectrum of the background is thermal, 
and that the excitation is not due to a local process, but rather a universal one. 
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Also included in this list for completeness are three late B and A 
type stars, that have recently been observed by other workers at the 120-inch 
Coude. Although T ^ has been left uncorrected for saturation m these cases, 
it is evident that the excitation here is also ~ 3°K. These three stars represent 
the first observation of CN in front of stars later than B . Since extensive H II 
regions are only associated with stars earlier than B 0, the very restricted 
ionization associated with stars of such comparatively late spectral type as these 
three provides very good evidence that collisions with charged particles do not 
materially contribute to T^. 

Hence the temperature invariance, while not conclusive, is to be 
considered to be very strong evidence that the observed interstellar CN excitation 
is due to a universal mechanism, which is an invariant property of the universe 
and independent of the CN location within the interstellar medium 

B. HI- CN Velocity Correlations 

For CN well embedded man HI region, we expect a velocity correla- 
tion between the CN absorption lines and the H I 21cm emission of this cloud. 

It should be noted that a HI- CN velocity correlation does not prove that the CN 
is in fact embedded in an HI region, since an H II region will also emit 21cm 
line radiation (Riegel 1967), however, it would seem'to be a necessary condition. 

Herbig (1968a) has shown that there is such a velocity correlation 
for the cloud in front of f Ophiuchi. We have also observed the CN and the H I 
21cm velocity profiles for Q Ophiuchi, as well as ihose for several other stars 
showing interstellar CN. The results are presented in Table 4, and Figures 
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25, '26, and 27 with the velocities referred to the solar rest frame. 

The H I 21cm radio observations were made with the 300 foot transit 
telescope of the National Radio Astronomy Observatory (Findlay 1963), -with the 
100 channel autocorrelation receiver (Shalloway 1964) at an angular resolution 
of 10'. 

The frequency resolution of the system’ for these’ scans was set 
between 2. 5 kHz and 5 kHz, depending upon the desired total frequency range 
for the scan. This corresponded to a velocity resolution of between 53 and 
1. 06 lcm/sec /channel. On the other hand the velocity resolution of the optical 
observations was ~5km/sec, and the optical absorption lines were barely re- 
solved. However, when we superpose the velocity profiles of CN, Ca II, and 
HI, 'we see a reasonable correlation. 

Figure 26 shows these profiles' in the direction of Q Ophiuchi. Note 
that 'the optical absorption lines have approximately the same velocity maxima, 
but that the’ 21cm ’counterpart to these occurs as an assymetrical bump on the 
side of the H I peak. This peculiarity was noticed by Her big (1968a). Presumably 
the H I emission originates both m front of and behind ^Ophiuchi, and represents 
two independent partially resolved clouds. That the 29km/sec Ca II cloud 
appears ‘to have no strong counterpart in either GN or H I is probably due to 
the rather high degree of saturation of the Ca II lines, not shared by the H I 
emission of the' CN absorption. In addition the highly irregular occurrence of 
CN absorption suggests that a lack of occurrence of CN absorption for this cloud 
should not be considered surprising. 



In Figure 27 the H I and CN velocity profiles are given for both 
Q Persei and X Persei, which are located within 1° of each other. Adams (1949) 
noted the identical interstellar velocities as well as other similarities m the 
spectra of these two stars. We now see that this similarity also holds for the 
21cm velocities as well. 

Finally Figure 28 shows the Ca II, CN, and H I profiles m the 
direction of 20 Aquilae. To the extent that the signal-to-noise ratio of the 
CN tracing permits, one can say that there appear to be three distinct clouds: 
the first at ~23km/sec containing Ca II, and possibly CN, a second at -12km/sec 
containing Ca II, H I, and CN, and a third at -2km/sec containing only H I 
and CN. 

The appearance of a CN feature at -2km/sec is somewhat curious, 
since no corresponding feature occurrs for Ca II. The fact that in this cloud 
H I is observed, but Ca n is not, might be explained through the assumption 
that this cloud is behind the star; however, to our knowledge there is no other 
case known of any interstellar cloud m which the CN absorption is stronger than 
that of Ca II, so that m this case the CN feature must only be noise. 

There is a good possibility that this CN feature is merely a dimple 
in the plate gram. On the other hand it is tempting to suggest that this is a 
unique cloud. It appears that only future observations with higher signal-to-noise 

ratio profiles for the CN absorption in this star can settle this question. 

These four stars are the only ones whose signal-to-noise ratio 

permitted a comparison of velocity profiles. The velocity maxima of these 
lines, as well as those of the CN lines in the spectra of several other stars, 
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are compared m Table 4 in the column labeled v. The entries m this column 
were made from traveling-microscope measurements of the CN lines on the 
spectral plates directly. 

All of the entries in this column (except possibly the aforementioned 
cloud in front of 20 Aquilae, whose existance appears dubious) shows a reason- 
able correlation of the molecular Ca II, and H I velocity. From these obser- 
vations (and from the previous work) we get the picture of discrete H I clouds, 
in which the molecules are well mixed, with the Ca II and H, as being the place 
of origin of the interstellar lines. 

Thus all of our data are consistent with the idea that the CN is well 
embedded in a normal H I region, and that the observed rotational excitation of 
the interstellar CN is not due to collisional or local mechanisms. 



CHAPTER 10 
* 

C ONCLUSIONS 
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In this dissertation T^(CN), the J = 0 -* 1 excitation temperature 
of interstellar CN, has been determined to varying degrees of accuracy in the 
direction of eight stars. The primary data for this analysis consisted largely 
of old plates found in the Mt. Wilson files, but included six spectra of Q Ophiuchi 
obtained by Herbig with the Lick 120-inch telescope and four spectra' of 
BD+66 0 1674 and BD+66° 1675 obtained’by Munch with the 200-inch telescope. 
This excitation temperature has also been determined in three other directions 
by other workers. A summary of these results is listed in Table 4 and discussed 
in Chapters 5 and 9. 

The important result of these measurements is that, 1 to within ob- 
servational uncertainties, T^(CN) appears to be independent of direction, as 
would be expected if the molecules were excited by cosmic background radiation 
alone. Unfortunately the uncertainty of T^(CN) is often of the order of 1. 0°K. 
New observational work on. the stars listed m Table 4 is therefore of high priority, 
and, of course, interstellar CN lines should be sought in the spectra of other 
stars. 

In the two cases with best signal-to-noise ratio, in which CN is ob- 
served against the 2nd magnitude stars Q Ophiuchi and Q Persei, T^ 0 (CN) is 
m remarkable agreement with the direct radio observations of the background at 
wavelengths longer than 3. 3mm (Figure 1). Consider in particular the series 
of direct measurements made by Wilkinson and his collaborators at Princeton, 
which represent the most extensive collection of direct observations made by a 
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single group, and those which extend to shortest wavelength. A weighted mean 

of their four measurements (extending from X 3. 2cm to X 3. 2mm) yields for 

the brightness temperature of the background radiation T = 2 68 ± 0. 11°K, 

j3 

while we find for COphiuchi T^(CN) = 2. 74 ± 0. 22°K, and for Q Persei 
T^q(CN) = 2. 74 ±0. 30° K. A graphical comparison of the radio observations 
and CN temperatures is shown in Figure 29. 

Although the invariance of T_^(CN), and the velocity correlation of 
CN and H I, are the strongest evidence that the CN lies in normal H I regions, 
a major part of this thesis has been devoted to considering the possible role of 
other interstellar processes in the molecular excitation. All of the evidence we 
have been able to assemble points to the conclusion that the CN molecules lie in 
normal H I regions where collisions or the scattering of photons are unable to 
account for 3°K excitation. 

We first considered radiative excitation in Chapter 6, and concluded 
that in general in the interstellar medium indirect radiative excitation mechanisms 
by way of fluorescent cycles are far too weak to account for the observed excitation. 

Next we considered collisional excitation, and assumed that the CN 
m front of Q Ophiuchi resides m the dense H I layer whose conditions are those 
found by Herbig (1968a). There is good evidence m favor of this assumption 
We find that ions may be neglected but electrons and H atoms may not. These 
latter two components provide roughly equal CN rotational excitation, however, 
using Herbig 1 s upper limit for the densities of these in this cloud, we find that they 
can contribute at most a small portion of the observed excitation temperature. 


This fractional upper limit is not as low as one might desire, however, 
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and it appears that for electrons at least, considerable improvement of this 
upper limit is possible by further observational work. In addition, it would be 
useful to know if the high H I density of this cloud is characteristic of the clouds 
associated with the other stars in our survey. In this respect observations and 
an analysis similar to that performed by Herbig for Q Ophiuchi should be exten- 
ded to these other stars. 

We next assume that the CN resides in an H II region. We find in 

_g 

this case that a modest density of only 1. 4 protons/cm ' (and electrons) is 
sufficient to account for the ~ 3°K observed CN excitation. However, we find 
that in this case the resultant excitation temperature is a sensitive function of 
the assumed proton density , and it is, now hard to account for the observed in- 
variance of T ^ m eleven stars. In addition three of these stars are of rather 
late spectral type, and have fairly limited H II regions associated with them. 

This fact casts strong doubt on the suggestion that the CN is in such a region. 

' V 

Thus on the basis of our current knowledge of interstellar conditions, 
the equilibrium between the CN molecules and ~3°K background radiation is 
altogether intelligible. However, the interstellar molecular spectra offer still 
more valuable information on the spectrum of the background radiation at 
wavelengths inaccesible to ground-based observation. The absence of interstellar 
lines from the CN J = 2 level,, and from the first excited states of CH and CH+ 
imposes low upper limits to the intensity of the radiation at X ~ 1 32, 0. 560 and 
0. 359mm. Since we are concerned only with upper limits to the background, other 
physical processes in the interstellar medium turn out to be of very slight con- 
sequence in the final analysis, and it is only necessary to consider the case of 
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C Ophiuchi where the lines can be observed with best signal-to-noise ratio. From 
the equivalent width measurements and upper limits listed in Table 4 and described 
in Chapter 5, we find that 

“ 15 ->2 “1 “1 

IJ2. 64mm) 3. 24 x 10 erg cm sec ster Hz 

Iv (l* 32mm) < 2. 75 x 10 1 erg cm 2 sec 1 ster 1 Hz 

1^(0. 56 0mm) <3. 91 x 10 14 erg cm 2 sec 1 ster 1 Hz 

-14 -2 -1 -1 

I (0. 359mm)<10. 0 x 10 erg cm sec ster Hz 

V 12 -2 -1 -1 TT 

t /n i r a v , o ,rT erg cm sec ster Hz 

I^{0. 150mm)<4 2 x 10 ° 

The relevance of these upper limits to the result of the only attempt 
so far to detect the background radiation above the atmosphere (Shivananden et al 
1968) was discussed in Chapter 5. It appears that their recent rocket observations 
of a large background flux with wavelength in the region of our upper limits must 
either be of radiation with a highly irregular spectrum, or of radiation which is 
quite local in origin, and does not extend as far as the -15Km/sec interstellar 
cloud in front of Q Ophiuchi. 

A final word should be said concerning some of the technical aspects 
of this dissertation. A large part of the data analysis consisted of devising a way 
of adding together a large number of spectra using a digital computer. These 
techniques are of potential value m a number of fields of astrophysics where 
signal-to-noise ratio is a principal problem. The techniques are described in 
Chapter 4. 

In Appendix 4 there is also a long discussion of the problem of optimal 
signal processing of digitized spectrogram tracings. Relations are derived for 
the optimal filter function, the best estimate of spectral line strength, and an 
error estimate for this line strength. 
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Finally, a side benifit of this considerable improvement ni signal- 

' ~ 1 ' „ 13 + , r ' i 

to-noise ratio, is the detection of C H with an abundance which is m agree- 
ment with the value found on the earth and Venus. 
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APPENDIX 1 


MATRIX ELEMENTS OP THE MOLECULAR HAMILTONIAN 

We wish to calculate the matrix elements of the molecular Hamiltonian 


K =B(K - A) + AAS 


(equation 3. 3) 


in the Hund's case (b) coupling scheme 

K = N + A, J =K + S 


(equation 3. 2) 


To calculate the matrix elements of the operator S^,, we must first transform 
it into the laboratory frame. 


S , = S' S(lq)(® _ VJ ) (aPY) (A1.1 

q = -1 q 

! 

S(lq) are the tensor components of the spin operator in the labor- 
atory frame, and 2D is the finite rotation operator transforming the laboratory 
frame to the molecular frame (Edmonds 1968, Chapters 4 and 5). The particular 
matrix elements of the rotation operator which appear in equation ( Al. 1) may be 
written as normalized spherical harmonics (Edmonds 1968 , equations 4. 1. 25 
and 4. 2. 7). 

-1 (1) (1)* ’ -q (1) 

(2D (ctSv) (apY) ~ (-1) sVq (apY) 

I \l/2 _ (A1 2 

S , may then be expressed as the scalar product of two commuting tensor oper- 
z 

ators 


-q 

s , 4f) S (-1) q S(lq) Y (W 

q=-l' 


= /"4tt/3 S(l) , Y(l) 

The matrix elements of the scalar product of two such operators 
may be written (Edmonds 1968, equation 7. 1. 6) 



Ill 


< AK'SJ |S ,lAKSJ> = 

Z 


(_1)K+ S + J jj g K , 

. IKS. 


<S !1 S(1) jlS) (AK' I 1 Y(l)|| AK) 

The first reduced matrix element in equation (Al. 4) has been evaluated by 
Edmonds (1957, equation 5. 4. 3), and the second is evaluated in the usual way 
by selecting the easiest to compute of the components (A K'MlY^lA K M), and 
dividing by the appropriate 3j symbol. Thus 

<AK'|lY(l)j|AK> = (AK'ojY lAK0>(-l) K '/K' 1 kY 1 

» \0 0 0 J ( 

The wave function . j AKM) of the rotor may also be written in terms 

(K) 

of the matrix elements S (a g y) of the finite rotation operator. Since 
l/o M\ ^ ^ 

Yi 0 (B-(3/4n) ’ (f3) equation (A 1. 4) then becomes 


<ak'||y( 1 )!Iak) =(-i) K ' + a /k' i kA 1 1 /3 

\0 0 0 / 8tt 2 \4tt 


3 \ 1/ ' 2 [(2K'+1) (2k+l) l 1 ' 2 


* TJT 


The integral over the three $ functions yields a simple product of 
two 3] symbols (Edmonds 1957, equation 4. 6. 2) and equation (Al. 5) sim plifies to 

(AK' |l Y(l)|jAK>(-l) K,+A (3/4rr) l/2 [(2K'+1) (2K+l)] l/2 /K' 1 K\ , A1 


^■A 0 A/ 

Combining these equations, we obtain the requisite matrix elements of the mole- 
cular Hamiltonian (equation 3. 3) in the Hund's case (b) representation. 

(A K'SJ 1 3C | AKSJ) =B[K(K+1) - A 2 ] 6__ ITr + AA(-1) MC ' +K+S+J 

K. IS* , , . 


x [ (2K'+1) (2K+1) (2S+1) (S+l) sf JSK 1 K' IK 

1 KS l- A 0 A 
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APPENDIX 2 

EVALUATION OF THE REDUCED MATRIX ELEMENTS OF THE DIPOLE 

MOMENT OPERATOR 

We wish to evaluate the reduced matrix element 
< AKSJa P C (1) |1 A'K'SJ'a' ) 

in the Hand's case (b) coupling scheme 
K = N + A , J = K + S. 


First we factor off the S dependence using Edmonds' (1968) equation (7. 1. 7) 

(AKSJa ||C (1) j! A'K'SJ'a) = 

(A 2. 1) 

(-1) K+S+J ' +1 [(2J'+1) (2J+l)l 1;/2 jK' J' s|(AKaP C (1) || A'K'd^ 


We evaluate, as before, by selecting one of the components, and 
dividing by the appropriate 3j symbol (Edmonds 1968, equation 5. 4. 1), whence 


(AKJaP C(l) I' A'K”J"a) = (-1) 


S+J’+l 


f(2J'+l) (2J+l)] l/2 jK' J' S\ ( AKaM = 0 | C n | A'K'ciM = 0) 
jj K 1 


(A 2. 2) 


K 1 K' 
0 0 0 


In order to calculate the matrix element in equation (A2. 2), we must 
first transform this component into a frame (primed) which is fixed with respect 
to the molecule. Using Edmonds (1968, equations 4. 1. 25 and 4. 6. 5) we get 


C 10 - E** («SY> O' 

q 


(A 2. 3) 


4D 


In a fashion similar to that of Appendix 1, since 0.^(8) = ®qo 


we can write 


< AKaM K =0 |C 1() 



q 


/EK'ciM^, =o) = (-1) A (3/4tt) 1 ^ 2 (1/8tt 2 ) [(2K'+1)(2K+1)] 
^qO ^JvO Sm P da ' d P' d ^’ < a 1 |a' > 


1/2 

(A2. 4) 
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The integral over the three 2D functions may again be evaluated by 
using ‘Edmonds' (1968,. equation 4. 6. 2). Integrating* and transforming from the 
Hund's case (b) representation to the representation m which the molecular 
Hamiltonian is diagonal, we have 


<AKja | C (X) II A'K'J'a'} = 7) u. 

Osf L * 


K'K 


U A 

K'K' KK 


(-D 


1 /p 

A+SfJ'+l [(2K'+1)(2K+1)(2J'+1)(2J+1)] ^ jK'- J' S 

/J K 1 


£ <a| C' | a' >/ K 1 K'\ 
q lq \-A q li ) 


(A 2. 5) 
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APPENDIX 3 

VIBRATIONAL FLUORESCENCE RATE 


The Einstein B coefficient for a vibration - rotation transition is 

given by 


B 


KSSid! [ j "Vi +tJ,,+i) Wf 


JV, J'v r 


3h 2 c 


(A3. 1) 


2J" + 1 


where r ° t ^, is the effective electric dipole moment for the transition 

v n J",v , J I 

( Herman and Wallis 1955). Since the P and R branches have approximately 
the same wavelength ,-we write the sum over rotational transitions as 


B 


v = 0, 1 


8tt 

3h 2 c 


:P 


vib-rot . 2 
v = 0, 1 1 


Using Oppenheimer's (1926) model, for an upper limit to B we consider the 
effective charge to be at the internuclear spacing, whence 


, (vib-rot) ,2 2 _ 

1 v n = 0, 1 1 M UJ 


_ Q Q 

— = 2. 4 x 10" Debye 


where jj. is the permanent electric dipole moment of the molecule, and B and 

e 

uu are molecular constants of its ground electronic state, found in Table 2 
© 

and Herzberg (1959, Table 39). Hence we calculate 

14 3 -1 -1 

B = 45. 3 x 10 cm Hz erg sec 

v - 0, 1 
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APPENDIX 4 

NOISE FILTERING AND STATISTICAL ANALYSIS OF ERRORS 

In this Appendix we consider the problem of the optimal measurement 
of the strengths of spectral absorption lines embedded in photographic gram noise. 
This is a very similar problem to one which has been considered m the field of 
communication theory, i. e. , the determination of the strength, position, and/or 
presence of radar echos in receiver noise. In the treatment which follows we 
shall make use of the literature of this field. * 

The basic problem is (1) to determine the most probable value of the 
strength of a spectral line formed in the photographic grain noise noise of a micro- 
photometer tracing of a spectrogram, and (2) to set error limits for this result 
Let us represent the intensity distribution recorded on the plate as a function of 
wavelength by 

y( X) = a( Xh BS( X) + n( X) (A 4. 1) 

where: 

a(XD is the continuum which is assumed to be known or measureable 
arbitrarily accurately, 

s(X) is a spectral line whose shape is also assumed to be known. 

(The validity of this assumption will be discussed later. ), 

B is the height of the absorption line to be determined (the minus 
sign signifies. absorption), 

n(X) is continuous gaus Siam (Goetz and Gould 1937, Fellgett 1951) 
random noise (assumed wide-sense stationary), whose auto- 
correlation function 

T 

R • (u) = / n*(X)n (u+X) dX (A4. 2) 

-T 

is assumed to be measurable arbitrarily accurately. 


* For a good review of this field see An Introduction to the Theory of Random 
Signals and Noise (Davenport and Root 1958), hereafter referred to as D &R. 
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Let us also represent by y'(X) our measurement of a spectrogram 
on a densitometer with a slit, whose transmission as a function of wavelength 
is given by h^.^X). Then* 


y'<« *[y* h slit \ < A4 - 3 > 

A similar and not unrelated problem which should also be considered 
is the removal of the slit scanning function h introduced by the densitometer. 
As long as the slit is narrow with respect to the most narrow spectral features, 
one can neglect the signal distortion introduced by its finite width. A 
common proceedure has been to use a wide slit function as a crude noise filter. 

This is not an optimal proceedure (for a rectangular slit function), however, and 
it will be shown that a statistical estimation of line strength is inherently tied to 
the problem of optimal noise filtering. The maximum likelihood estimate of a 
line strength will be shown lobe a function of certain parameters of an optimally 
filtered signal**. Our analysis will treat first the problem of processing the 
direct signal y(X) , and then return to the problem of including the effect of a 
finite slit width. 


A. Maximum Likelihood Estimation 

Asa starting point we shall use a maximum likelihood estimation 
for the most probable value of (3 (VonMises 1964, p. 547), and set our error 
limits according to the criterion of the Ney man- Pear son theory of hypothesis 
testing (VonMises 1964, p. 504). 


* Here a*b represents the convolution product given by 

[ a*b] = P a(u - X) b(X)dX 

** Fellgett (1953) has discussed the filtering of densitometer tracings of spectro- 
grams, and come to the conclusion that the best filter is a "matched" filter. 
However, his filter is only derivable for an infinite length spectrogram, and 
he provides no suggestion for an estimate of spectral line strength. 
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If y =' ! y. y o ■ , y is a given set of measurements*, from which 

we would like to estimate 'the parameter p, then the aposteriori probability density 
that our measurement will have occurred is 

q(y) = p(y 1 P) p( P ) (A4. 4) 

where p(3) is the apriori probability density for the occurrence of j3. 

The apriori probability density p((3) presents a problem because it 
is unknown. To find it, one would need a large ensemble of identical' experimental 
situations, however we do not have such an ensemble. We do know, however, 
that 3 s 0, (i. e. we have no emission lines). Thus the apriori probability of an 
occurrence of 3 which does not satisfy this inequality is zero. 

j * 

Outside of this restriction, p(3) is unknown, and the usual practice** 
originally due to Gauss (1823) is to assume p(3) to be constant m this region. 
Equivalently 3 may be regarded as an unknown positive "constant". 

We define 3 to be the maximum likelihood estimate o f 3 equal to 
the value of 3which maximizes p(y|3) (called the likelihood function) . This is 
just equivalent to saying that we choose the value for 3 which makes our obser- 
vation the most probable one. *** It is generally convenient to maximize 
log P(y | 3 ) instead of p(y | 3 ). Since the logarithm is a real increasing function 
for positive arguements, the proceedure is equivalent to that of maximizing 

p(y 1 0)- 

For an error estimate we must define an interval H, which for a 
given probability Q , the actual value of 3 lies within this interval. Usually 
m physics when gaussian statistics apply, the domain of a parameter is infinite 
or very great m comparison with the error limit range. Then the most probable 

* For example the y^ may be the digitized intensity points of a spectrogram, or 
as we shall use later, the coefficients of an orthoganal expansion of y. 

** Further justifications for the proceedure is given by VonMLses (1964, p. 496, 
548). 

*** P(y I 3) is tacity assumed to possess one and only one maximum in the region 
where p( 3) 4 0. 
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value lies at the center of the interval and one says, for example, that for 95. 5% 
probability the true value of p lies in the interval between p - 2oand P +2a [or 

A 

with a probability Q = erf (n/j2 ) , the true value of plies in the region p- noand 

A A 

P + no] , and the most probable value of p = p . We shall adopt the 95. 5% 

A 

probability (n = 2), but find that p<na', the region is no longer symmetrical 
about P 

The probability that the true value of p lies within the region H is just 


where: 


Q(H) 


/ H P(y«=A| p ) P(p) dp 
4+H p(y cA lP)pU')^ 


(A4. 5) 


H is our error limit region 

A is the region of y for which we will decide that p H and, 

H+H is the allowed domain of p(p>0). 

Since p(P) is constant over H+H equation (A 4. 5) reduces to 
/ p(yc A | P) dp. 

Q (H ) = (A 4. 6) 

^H+H P(y ^ A i g)dP 


What we wish to find is the minimum size region for H such that for a given 
probability Q the actual value of p is in H. 

If n(A) is a gaussian random process, it can be shown (use is made 
of the results on pages 336 and 342 of D & R) by a rather involved analysis which 

4 

will be omitted here for the sake of brevity (and to keep from having to type all 
of the equations) that p ( X) p(y I p ) is of the form 


(.exp . [-(P- p ) 2 / 2 a 2 
p(P)p(y|p)=N m p 


] for p> 0 
otherwise 


(A 4. 7) 


where N is a normalizing factor. 
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f 




A A 

L ry(» -a(!0 ]£(i) di 

V* 


V A 

/. A s(X)((X.)dX 


(A 4. 8) 


a 8 A^+A 


L A S(X) f u> ^ 

and f{ X) is defined by the following integral equation 


(A4. 9) 


A +A 

J R (u-X.) f( X) dX = s(u) 
JL-A n 


(A4. 10) 


Now consider the troublesome domain of (3. To maximize the integral 
(A4. 6) and to find the maximum point of p(y | 3), we must consider the three 
cases which are illustrated in Figure 30. We must find an error limit region 
H of minimum size such that, with probability Q = erf (n// 2 ), 3 will be within 
H (our error criterion will yield the probabilities corresponding to the standard 
± ncr criterion used in physics). 


B. Case A 


If 3 s ncr* = <5 A erf 1 f Q/ 2-Qi 
m 


(A 4. 11) 


then 3 = 3 m and it is easy to show that a minimum error probability occurs for 
limits symmetrical about 3- Then with probability Q the true value of 3 lies in 

A *• A 

the interval 3 ~ P - 3 ^ 8+ t|. 

As can be seen from Figure 30, in this case we have for equation 

(A 4. 6) 

3 +p 

m 

L P(y I 3) d 3 

m 


/ p(y i 3) dp 
o 


erf r\//2 ~Q = 


(A4. 12) 
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Equation (A 4. 12) may be solved for rj to yield' 

n= a/2 erf 1 { T 1+erf ( P m /a/2 ) TQ/2^ . (A4. 13) 

We can also see now the reason for the definition of r|d. It is 

defined so that the lower limit of the integral in the numerator of (A 4. 12) is 

greater than or equal to zero i. e. - T) s 0. 

The limiting case occurs at (3 = ri, and from (A4. 13) we may solve 

-1 m 

for ri ; thus rp' = r[ = p =c/2 erf [Q/2-Q1. 

m 

For n = 2, 2 o' = 1. 7116a. 

Note that r\ reduces to the usual case encountered in physics for 

B well into the domain of 8, then 
m 

lim rj — a 

S3 -* a» 
m 

C. Case B 

If 0 < 8 ^ d , then 8=6 is still the most probable value of 

m m 

8 , but it is no longer in the center of the error region, and equation (A4. 6) 
becomes 

B + n 

J 0 P(y | P) d8 

erf n// 2 = Q = — — (A4. 14) 

I P(y 1 8) dp 

0 

Solving (A4. 14) for r| we get 

r] = a/2 erf 1 [ Q -(1 - Q) erf ( j p j /a/2) ] (A4. 15) 


D. Case C 

If P^ £ 0, the maximum value of p(y | p) occurs for p= P = 0, 
since p(y | p) is never stationary m the allowed domain of Pas can be seen 
from Figure 30. Then the lower limit of the region is zero, and the upper limit 
is again given by equation (A 4. 15), as long as we include the absolute value 
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sign, and define n as is shown in Figure 30. 


E. Interpretation 

These results can be interpreted simply m terms of optimal linear 
"matched 1 ' filters ( D &Rp. 244, Zadeh and Ragazzini 1952, and Van Vleck and 
Middleton 1946). 


. Let us define 

g(A L - X) = ( p (X) for A l ~ A ^ X ^ A^ + A 
1 0 otherwise 


(A 4. 16) 


CO 

y"(X) = - / T y(u) - a(u)1 g(X- u)du = T (y-a) * gl (A4. 17) 

and 


s"(X) = J s(u) g(X-u)du = [s*gl (A 4. 18) 

A 

Then 8 and a o are just given by 
m p 

= / s,t (A4. 19) 

u 2 = l/s" (A ) (A 4. 20) 

p -Lt 

It is clear that y"(X) is a representation of y which has been filtered 
by the function g. If the value = y"(A^) s 0, it is just the maximum 

likelihood estimate of the signal strength. This is equal to the height of the 
filtered, zero restored, signal at the line center, if a normalized filter is used. 

The integral equation for the filler function may be written by com- 
bining (A 4. 15) and (A 4. 16), whence we get 


s^- u) = / 


R (X- u) g(X) dx- 

a l -a 


for -Aj£u ^A 


= JH n *g \ = s(u) 


(A4. 21) 
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Here we have, defined 

s’ (u) = s (A^ - u) 


F. Rem o val of Effect of Finite Slit W i dth 

Let us now return to the problem of inverting the slit function. In- 
stead of using the function g to filter our data, we will use the function h , 
whose purpose is to both remove the affect of h , and to apply g. The tandem 
filter function for two linear filters applied sequencially is given by the convolution 
of the two individual filter functions (Cheng 1961, p. 245), 


h * h = g 
c slit 


(A4. 22) 


Convolving both sides of equation (A4. 22) with R* h ^ , where 
h ($ — h '(-\)» and using equation (A4. 21) we have 

Silt Silt 

(R * h *h . J *hc =R *g *h 
n slit , slit n 


slit 


(A 4. 23) 


= s’ h = s' 
slit 


where we have defined 


s'U> “ r 8 * h s ii t \ = s '(A l -A) 


(A 4. 24) 

By the use of the Wiener -Khintchin theorem ( D & R p 183) we can 


identify 


R * h * h ... =R» 
n slit slit n 


(A4. 25) 


as the autocorrelation for noise which has been filtered by h ^ ; similarly s' 
is just the pure signal filtered by h . Thus we have 


[R’ * h 1 = s'( A - u) 

n c~u 


(A 4. 26) 


This will just be noted as the integral equation for the optimal filter for extracting 
a signal, which has been modified by convolving the slit function with it, from noise 
which has been modified in the same way. Thus all we need for the calculation of 
the optimal slit and noise removing filter function is R* , not R This is 


n 


n 


fortuitous in that R' may be measured directly from the digitized tracings by 
n 
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a simple numerical integration, where as JEt would require additional .optical 

n 

measurements. 


G. Error Estimates for Spectral Line Depth 

A 

From equations (A4. 22) and (A4. 24) an error estimate for (3 
(taking into account the three cases previously described) is given by 


ct = [h * h * R \ - B*(0) 
p n n n 0 n 


(A 4. 27) 


where 


\ M - 




(A4. 28) 


c 

is the normalized filter which applied to y yields y n (normalized), then R' is 

■* n 

just the autocorrelation function of pure plate grain noise which has been acted 

2 

upon by h . Thus a measurement of a 0 becomes very simple: just measure 
n p 

R^(0) i. e. the autocorrelation of filtered noise in an adjacent region of the plate 
where no signal is present. This is the way in which errors were calculated for 
Table 4. 



124 


APPENDIX 5 

CALCULATION OF FILTER FUNCTIONS 

Integral equation (A4. 26) becomes particularly easy to solve if one 

assumes that the noise spectrum is a constant (white noise), for which R^(X ) 

= W 6(>) and W is given by 
n n 

00 

W = / R' (X) dX (A5. 1) 

n J n 


Then we have the immediate solution 


g(X) 


s^- X) 
0 


for /^-A^X^Aj+A 
otherwise. 


(A5. 2) 


and our likelihood test becomes simply a cross-correlation of the measured 
signal with the assumed signal. 

In practice, RI (X ) will be highly peaked at the origin with a peak 
width controlled by the slit width. The measured peak half widths are typically 
. 015A. 

Using autocorrelation functions R^ of similar form to those measured, 
equation (A4. 26) has been solved for a gaussian signal shape with a line width 
parameter . 05 A, with a method described below. The equation for an optimal 
filter function is of the form 


where 


A^ 

/ R' (u - X) f(X) dX = s(u) 
-A 


R' 

n 


R' 



(A5. 3) 


.In the case when R' (X.) « 6 (X)- we know immediately a good approx- 

\ E- 

imate solution is 


h(X) « s(X) 


(A5.4) 


Denote the linear integral operator of (A5. 3) by £(h‘(X) ) then (A5. 3) becomes 
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£(h(X) ) = s(X) 


(A5. 5) 


let h = s + where ^ «h, then £(h ) = s - £(s) = F(s) = s ± and s 1 «s. 

We now have a new integral equation to solve for the correction: 

£(h x ) = s i (A5.6) 

By iterating the same proceedure, we get for a solution 


where 


00 



s^ = s — £(s) = F(s) 

s 2 = s ^ — £( =F(s 1 ) = F 2 (S) 


(A5. 7) 


(A5. 8) 


s 

n 


n-1 


-£(s J 

n-1 


F(S n-l ) = fU(s) 


For the cases calculated, three iterations were sufficient to reduce the correction 
to less than . 1%. 

Filters were calculated using this technique assuming that the shape 
of a spectral line is gaussian (Stromgren 1948). Although this assumption has 
been criticized by Munch (1957) and Blaauw (1952), the actual velocity distribution 

j 

within an interstellar cloud, and hence the resulting spectral line shape is 
peculiar to the cloud under consideration. Since most of the other high signal- 
to-noise ratio lines in the stars of our survey are well fit by a gaussian line 
shape, this assumption appears to be very good. 

T-he difference between the assumed gaussian signal and the calculated 
filter function is found to be only about 5% (actually the solution closely approximates 
a gaussian slightly narrower). Since the signal shape is not known this accurately, 
a gaussian filter function was used for the processing throughout. 
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TABLE 1 

IMPORTANT RESONANCE TRANSITIONS OF INTERSTELLAR MOLECULES 


Optical Transitions 


Molecule Electronic 

Lines Origin- 

Excited State 

X air (A 0 ] 

Ground 

ating in that 

levels 


State levels 




CN (X 2 ^) 


(B 2 S + ) 


J = 0 

R(0) 

J = 1 

3874. 608 

W 

II 

R(l) 

J = 2 

3873. 998 


P(D 

J = 0 

3875. 764 

II 

R(2) 

J = 3 

3873. '369 


P(2) 

J = 1 

3876. 312 

CH <X 2 II) 


(A 2 A) 


KJ=1 1/2 

Rzd) 

KJ=2 3/2 

4300. 321 1 

KJ=1 3/2 

Ri(i) 

KJ=2 5/2 

4303. 948 


R Qi2(i) 

KJ=2 3/2 

4303. 621 


S R 2 i(i) 

KJ=3 5/2 

4287. 620 

KJ-2 3/2 

Rj(2) 

KJ^3 5/2 • 

4296. 636 


Q^2) 

KJ=2 3/2 - 

4312. 705 


Q Rl2<2) 

KJ=25/2 

4313. 027 

CH + (X J S + ) 


(A *11) 


it 

o 

R(0) 

J = 1 

4232. 535 C 

J = 1 

R(D 

J = 2 

4229. 336 


Q(l) 

J = 1 

4237 555 




Microwave Transitions 

Observed Intensity Ratio 

s 

Allowed Transitions X 

within ground state (mm) 


R(1)/R(0) 

P(1)/R(0) 


*2 

1 

J=0 - J=1 

2. 644 

R(2)/R(l) 

P(2)/R(l) 


*3/2 

1 

J=1 - J=2 

1. 322 

Rl(l)/R2d) 

(9/5) d 

*1. 525 

KJ=1 1/2 - KJ=1 3/2 

0. 560 

R Q I2 (1)/R2(1) 

(1/5) 

.392 



S R a (l)/R2(l) 

( 0 ) 

. 083 



Rg(2)/Ri(l) 

Qa(2)/Ri(l) 

(20/27) 

(1/3) 

* . 820 
. 267 

KJ=1 1/2 - KJ=2 3/2 

0. 150 

Q R 12 (2)/R 1 (1) 

(1/27) 

. 0224 



R( 1)/R(0) 
Q(1)/R(0) 


* 3/2 
3/2 

J=0 - J=1 

0. 359 


Jenkins and Wooldridge (1938) 
Gero (1941) 

Douglas and Herzberg (1942) 
Hund's case (b) 
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TABLE 2 

MOLECULAR CONSTANTS 



CN 

CH 

CH+ 

-1 

a 


c 

B (cm ) 
-I 

1. 8909 

14. 457 
„ _ c 

13. 9318 

A (cm ) 
M D ) 

** 

1. 45 ±. 08 

27. 95 

d 

1. 46 ± 0. 06 



a Poletto and Rigutti (1965) 

k Thomson and Dalby (1968); note that,^ = 1. ID (Arpigney, 1964) was used for 
the calcuations of Chapter 8. 

e Herzberg (1950) Table 39 

d Phelps and Dalby (1966) 
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TABLE 3 

DATA ON PLATES USED IN NUMERICAL SYNTHESIS 


Star 

Plate No 

Date GCT(beg ) 

Exposure 

(minutes) 

Slit 

Width 

Emulsion 

Camera 

(Grating) 

Spectral Seeing 
Width 

Observer Fraction of Synthesis 
i CN CH CH* 

f Oph 

Ce2313 

5/13/40 10 55 

240 

6 

E-33 

114" 

1 0 

2-3 

A 

0 

0 

3 1 

4 0 


Ce2317 

5/14/40 10.31 

270 

6 

Cr Cont 

114" 

0 5‘ 

3-2 

A 

0 

0 

4 1 

16 0 


Ce2328 

6/15/40 7 39 

360 

6 

Cr Cont, 

114" 

0 7 

1-2-3 

A 

5 

5 

0 0 

o a 


Ce2334 

6/17/40 7 45 

245 

5 l/Z "*■ 

Cr Cont 

114" 

O 5 

3-4 

A 

$ 

0 

0 0 

0 0 


Ce2381 

7/26/40 9: 08 

127 

6 

Cr Cont 

114" 

0 8 

3-2 

D 

28 

9 

23 6 

I 0 


Ce2414 

9/13/40 8 03 

50 

7' 

Cr Cont 

114" 

0 5 

2 

D 

0 

0 

10 3 

1 9 


Ce2495a 

1/17/41 16 48 

’ 18 

7 * ' 

Cr Cont 

114" 

0 2 

3 

D 

0 

5 

0 0 

0 1 


Ce2495b 

1/17/41 17 06 

44 

7 

Cr Cont 

114" 

0 2 

3 

D 

1 

7 

0 0 

1 4 


Ce2516b 

3/16/41 15’ 19 

10 

7 

Cr Cont 

114" 

0 2 

3 

D 

0 

0 

0 2 

0 0 


Ce2516e 

3/16/41 15 29 

48 

7 

Cr Cont ' 

114" 

0 2 

3 

D 

7 

3 

1 7 

4 9 


Ce2522 

3/17/41 13 04 

63 

7 

Cr Cont 

114" 

0 3 

3 

D 

1 

4 

3 0 

0 0 


Ce2538 

4/16/41 11 42 

106 

8 

Cr S C 

114" 

0 3 

4-1 

D 

0 

0 

4 6 

3 6 


Ce2542 

4/21/41 10 09 

365 

6 

Baked 

114" 

0 7 

3-1 

A 

5 

2 

12. 7 

0 0 






Process 










Ce25S9b 

5/11/41 11- 26 

30 

8 

Exp 144-0 

114" ' 

0.2 

2 

D 

3 

3 

0 0 

0 0 






164,503 










Ce2573 

5/27/41 9 59 

342 

6 

Cr. Cont 

114” 

0 7 

2 

A 

0 

6 

4 7 

2.3 


Ce2578 

6/ 4/41 8: 16 

360 

6 

Cr. Cont. 

114" 

0 5 

4-2 

A 

3 

1 

4 4 

2.9 


Ce2581 

6/ 5/41 8 50 

332 

61/2 

Cr Cont. 

114" 

0 5 

4-1 

A 

3 

1 

0 0 

0 0 


Ce2582 

7/ 4/41 7 39 

231 

6 

Cr Cont 

114" 

0 3 

6-5 

A 

4 

1 

14 7 

21 2 


Ce3847 

6/24/45 7: 45 

35 

91 

103-0 

114" 

0 8 

2-3 

A 

1 

6 

0 0 

2 B 


EC28, 

3/15/60 




160" 



H 

4 

4 

0 0 

o a 


EC28; 

3/15/60 




160" 



K 

3 

2 

0 0 

0 0 


EC38j 

3/18/60 




160" 



H 

0 

0 

4 1 

1 7 


EC38 2 

3/18/60 




160" 



H 

0 

0 

6 1 

18 3 


EC55, 

4/16/60 




160" 



H 

0 

9 

0 0 

0 0 


EC55j 

4/16/60 




160" 



H 

1 

5 

0 0 

0 0 


EC220 

8/ 8/60 




160" 



H 

2 

1 

0 0 

0 0 


EC162, 

7/ 4/60 




160" 



H 

4 

1 

0 0 

5 4 


EC162 2 

7/ 4/60 




160" 



H 

3 

7 

0 0 

4 1 


ECI890 





160" 



H 

0 

3 

0.0 

0 0 


EC1969 

6/12/63 




160" 



H 

5 

5 

3 1 

8 4 


? Per Ce2432 

10/13/40 

10 57 

50 

96 5 

Cr H S 

114" 

0 7 

4 

A 

5 3 

Ce2480a 

12/20/40 

7-15 

40 

7 

Cr Cont 

114" 

0.2 

3 

D 

17.3 

C e24 80b 

12/20/40 

7 55 

30 

7 

Cr Cont 

114" 

0 2 

3 

D 

10 1 

Ce2480c 

12/20/40 

8 25 • 

20 

7 

Cr Cont 

114" 

0.2 

3 

O 

1 5 

Ce2481 

'12/20/40 

9 33 

150 

7 

3J + Am 

114" 

0 2 

2 

D 

3 0 

Ce2492a 

1/17/41 

8 18 

36 

7 

Cr Cont 

114" 

0.3 

4 

D 

41 2 

Ce2492b 

1/17/41 

8 54 

20 

7 

Cr. Cont. 

114" 

0 3 

4 

D 

21 4 


55 Cyg 

Ce2378 

7/24/40 

11 20 

290 

6 

E-33 

114" 

0.5 

3 

D 

59 

0 


Ce2390 

8/18/40 

10 07 

378 

9 

Cr. Cont. 

73" 


3 

D 

30 

0 


Cel3793 

8/18/60 





73" 



H E B 

11 

0 

AE-Aur 

Cel261 

12/ 6/36 

9:22 

150 

91. 5 

m-F+Trl 

32" 


4-3 

M 

33 

0 


Cc3658 

12/ 6/44 

10. 04 

110 

91,5 

103-CBakd 114" 


5 

A 

10 

0 


CS9719 

3/ 5/55 

7 33 

51 

114 

IIa-0 

32" 


2-3 

SI 

57 

0 

20AqI 

Ce3554 

9/ 4/44 

8 21 

90 

92 

103-0 

114" 

0 5 

3-4 

A 

100 

0 


X Per 

Cel534 

10/22/37 

13 

40 

90 

12 

E-33 

32" 

5 

Sa 

20 0 


Ce2867 

9/24/42 






114" 

0 3 

A 

80 0 

BEH- 

66°1675 

PC6056 

8/ 1/61 

13 

03 

167 

2/40 

Ha-OBaked 

36" 

2-3 

Mu 

78 0 


PC6131 

8/31/61 

12 

12 

180 

4/23 

lla-OBaked 

36" 

3 

Mu 

22. 0 


BD+ 

66°1674 

PC7520 

9/ 5/63 9 32 200 

2/50 

Ila-oBaked 

36" 

2-3 Mu 

46 0 


PC909 1 

12/263/65 6 156 3596 

2/45 

lla-OBaked' 

36" 

4 Mu 

54 0 

Ce = Mt 

Wilson 

6 05 319 

100-mch, PC = Mt. Palomar 

200-inch, 

EC = Lick 

120- inch 

& 1-2 
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TABLE -4 


bl ir 

'} 

MK 

b" 

C) 

i" 

T) 


b 

Mol/ 

Line 

Equivalent 

i i . 

r 

Al" 

T(* K)' ,k 

vlKm/stc) 1 



1 vpc 

{Km/scc) Aiom 


Widlh(mA) 1 

* 

('K) 

t 0|)h 

2 »b 

09 )V 

+23 b~ 

S~3 



HI , 

21cm 





-12 7141 • 






3 

54® 

CAII 

K 





-29 0 e I 0 3/6] 







61® 

CN 

R(0) 

6 621, 29 

572 ± 032 



-14 8^1 0 4/6) 








CN 

R(D 

2 06 £ 29 

155 ± 024 

20 

2 72± 22 

-13 8^[ 1 0/4] 








CN 

P(l) 

1 12 L 29 

082 ± 022 

25 

2 81 ± 31 

-14 9 e [l 1/31 








CN . 

R(2> 

(0) < 44 

(0) < 032 

10 

(0) <5 49 









CH 

Rj(i) 

13 40+ 03 

1 50 ± 03 



-14 6«[0 45] 








CH 

Ri(l) 

( 11} < 53 

< 008) < 039 

71 

(4 54)<C6 33 









CH 

R^?) 

( 11) < 53 

( 008) <039 

1 8 

<19 4)<28 5 









CH + 

R{0) 

19 60 £.08 




-14 0®[0 7/6) 








CH + 

R(l) 

( 14) < 72 

< 009)0 048 

1 07 

<6 86) <9 01 






* 


. 

C 1S H* 

R(0) ‘ 

60 ± 31 

044 + 023 
- 013 

c .2/ c l 5 = 62 + _ 188 


£ Per 

2 83 

Bllb 

-IS 7 

162 3 


61 

HI 

21cm 





13 4[2t 








CN 

R(0) 

9 12£57 

079 
859 - 070 



12 31 1 3/2] 








CN 

R(l>' 

3 02 £ 57 

234 ± 049 

30 

2 72 £ 30 









CN 

P<1) 

1 24 £57 

092 ±044 

30 

2. 43 ± 54 

• 








CN 

R<2) 

{ 44) < 1 23 

{ 032)<091 

29 

<4 55)<8 05 









CH6CH+ 

, 

« 


13 8 f [7l 








CAII 

H&K 





13 5 f 

S3 Cyg 

4 83 


B3Ia 

85 7 

1 

0 

CN 

R(0) 

b 6 i 3, 0 

32 £ 16 











CN 

H(l) 

( 2 0S> C 4 05 

( 09) < 23 

14 

(2 8)<5 5 


AEAur 

5 30 

08 5v 

- 2 3 

172. 1 

1. 

0 

HI . 

■21cm 





2 5, 11 5(broad) 








OH6.CH+ 



















15 7)[ 6] 
14 5 f 








CAII 

H&K 












cn 

R(0) 

5 90*2.90 

28 £ 15 











CN 

R(D 

2. 67 *2. 57 

12 £. 11 

. 2 

3 5 ±2, 3 


20 Aql 

5 27 

B3Iv 

- 8,3 

28 2 

1 

0 

HI 

2l<-m 





-11 5, -2 0 
-22 7 f , -12 l f 








CA 

K 












CHLCH+ 





-13 4 r | 3] 









- 












CN 

R(0) - 

22 9 ±8. 0 

1 64 +1 « 
- .78 



-12. 2,-2. 8fl 2/3] 











+ .53 


+1 8 









CN 

R(i) 

7 6 ±7. 2 

38 - 35 

5 

25 -15 



HD 

11953 

5.68 Alla 

- 3 0 

132 9 

CN 

R(l> 




3 7£ 7° 


HD 

208501 

5 79 B8Ib 

1 7 

100.4 

CN 

R(l) 




2 8£ 4 d 


X Per 

6 08 

-17 1 

163 1 10 

HI 

21cm 





13 3, 22 8(1] 
15 4 f 





CAII 

H&K 









CN 

R(0) 

18 4 ±3 1' 

+ 33 
1 15 - 26 



13 5, 23 2(1] 





CN 

R(l) 

6 8 ±3 1 

34 t 15 

- 4 

2 8± 8 






CH&CH+ 





-15 2 f l6] 

HD 

BS 










2057 1 

6* 10 IJ-III 

-20 6 

172 4 

CN 

R(l) 




~3 a 


BD+ 

k 










66* 1675 

0. 05 07 

+ 48 

118 0 5 

HI 

21cm 





-14 2,-19 5[ 1] 




• 

CAII , 

K 





-15 iSf 2/2 1 





CH 

Rid) 





-21 3?( 7/2] 





CN 

R(0) 

125 5 £14 0 

1 95 ± 47 



-17 2S(1 0/1 SI 





CN 

R(l) 

39 6 ±14 0 

40 £ 16 

56 

2 39 £ 43 


BD+ . 
66* 1674 

9 5 0 

+ 48 

118 0 5 

HI 

21cm 





-14 2,-19.511) 





CAH 

K 





-17 9S[ 7/2] 





CH 

Rzd) 





-17 6B[1] 





CN 

R(0) 

55 9+7 3 

60 £ 10 



-17 4^(1] 





CN 

R(D 

13 8 ±7. 3 

. 13 £ 07 

. 15 

2 45 ±63 




NOTES TO TABLE 4 


G. Herbig, Lick 120-inch Coude, private communcation 

From spectra taken by G. Munch (1964) with 200-inch Coude 

« 

M. Peimbert, 120-mch Coude, private communication 

V. Bortolot Jr; and P. Thaddeus, 120- inch Coude, private communication 

Measured by G. Herbig (1968a) 

Measured by W. S. Adams (1949) reduced to the wavelength standards used 
by Herbig (1968a) , Munch (1957) 

Measured by G. Munch (1964) 

i 1 

A T is the correction (decrease) for line saturation 
» f 

Equivalenth width, optical depth, and temperature error limit ranges corres- 
pond to 95% confidence level except when aline has not specifically been detected 
with this confidence level. Then the upper limit given is to 99. 7% and the 
number in parenthesis is the most probable value. 

Velocity errors follow the notation of Herbig (1968a) with the first number in 
brackets being the standard error and the second number being the number of 
plates measured. 

For a given pair of levels i and j the temperature calculated from n /n. is 
listed beside optical depth t.., . n. is always calculated from the strongest 
line originating m level j. 
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Figure 1. Direct measurements of the brightness temperature T of the 

B 

fireball radiation. These are numbered in chronological order. 



Wa.|ele^gth 

t b ( °K) 

Investigators 

(11) 

73. 5 

3. 7 ±1. 2 

Howell and Shakeshaft (1967) 

(12) 

49. 2 

3. 7 db 1. 2 

Howell and Shakeshaft (1967) 

( 4) 

21 0 

3. 2 ± 1. 0 

Penzias and Wilson (1967) 

( 3) 

20. 7 

2. 8 ± 0. 6 

Howell and Shakeshaft (1966) 

( 1) 

7. 35 

3. 5 ± 1. 0 

Penzias and Wilson (1965) 

( 2) 

3. 2 

3. 0 ± 0. 5 

Roll and Wilkinson (1966) 

( 1) 

3. 2 

2 69 + °- 16 
^ 59 -0. 21 

Stokes, Partridge and 
Wilkinson (1967) 

( 8) 

1. 58 

' 2 78 +°* 12 
2 - 78 -0. 17 

Stokes, Partridge and 
Wilkinson (1967) 

( 5) 

1. 5 

2. 0 ± 0. 8 

Welch, Keachie, Thornton, 
and Wrixton (1967) 

( 9) 

0. 924 

3. 16 ± 0. 26 

Ewing, Burke and Staelin (1967) 

( 6) 

0. 856 

2- *6 + H 
- 0. 22 

Wilkinson (1967) 

(10) 

0. 82 

2. 9 ±0. 7 

Puzanov, Salomonovich and 
Stankevich (1967) 

(13) 

0. 33 

2.46 + ^ 
- 0. 44 

Boynton, Stokes and 


Wilkinson (1968) 
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Figure 2. Term diagram which is required for an atom or molecule to be able 
to be used to measure the cosmic microwave background with the techniques 


described here. 
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Figure 3. Partial term diagrams for molecules observed optically in the inter- 
stellar medium showing transitions which produce the strongest interstellar 
lines - All levels shown, m both the ground and excited electronic states, are 


in the ground vibrational state. 



f) 


Rote) 


R(0) R(! ) 
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Figure 4 . Schematic diagram of equipment used for tracing plates at California 
Institute of Technology. 




MAGNETIC TAPE 

ROBINSON LABORATORY WILLIS BOOTH COMPUTING 

OF ASTROPHYSICS CENTER 
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Figure 5. Schematic diagram of the areas traced on each plate showing placement 
of razor blade for fiducial mark. Arrows show direction of scan and slit orien- 


tation. 



TO PHOTOMULTIPLIER 




STELLAR 

SPECTRUM 


COMPARISON 

SPECTRA 


CALIBRATION 

WEDGES 


0 — ► INDICATES DIRECTION OF SCAN 
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Figure 6 Wedge bars and corresponding intensities resulting from tracing 
plate perpendicular to stellar spectrum. Note that opposite sides of the wedge 


were put on at different intensities. 



TRANSMISSION 
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Figure 7. Wedge after compensation for uneven intensities on opposite sides of 
plate. Square shows average continuum intensity. 




INTENSITY 



Figure 8. Plot of wavelength versus digitization index with least square 


straight line fits. Different initial wavelength resulted from fact that razor 
could not be placed exactly perpendicular to the stellar spectrum. Numbers 
refer to symbolic identifications of iron lines. 



I LI 


3880 


3878 


3876 


t 3874 


3872 


3870 
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Figure 9 . Syntheses of the spectra of several stars showing absorption of 
interstellar CN. All but £ Gphiuchi and £ Persei have been numerically 
filtered. Vertical magnification factor and number of spectra for each case 
are noted at left. Stellar lines are evident in the spectra of several of these, 


notably Q Persei and 55 Cygni. 
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Figure 10. Synthesis of 17 spectra of £ Ophiuchi showing absorption by inter- 
stellar CH. Spectrum has heen filtered by a gaussian filter function of width 


. 05 A'. 








Figure 11. The result of adding together 19 spectra of £ Ophiuchi showing 

absorption by interstellar CH+. Spectrum has been filtered by a gaussian 

o 13 + 

width of . 05A. The location of the isotopically shifted R(0) line of C H 


is also indicated. 






Figure 12. Limits on the intensity of radiation in the interstellar medium set 

by this work (M W) and by further observations by Bortolot et al (1969 - L). 

-9 -2 -1 

The area under the line AA ’ represents the flux 5x10 W cm ster 
observed by Shivananden et al (1968). 
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400 
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600 
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2 2 

Figure 13. Typical fluorescent cycles through Iland 2 
level J with respect to level J". Note for cycle through 
(Q branch) transitions are allowed. 


states which populate 
II state, AJ = 0 




Figure 14. Electronic term diagram for CN showing the well known red and 


violet systems and the ultraviolet systems studies by Douglas and Routly 
(1954, 1955) and Carroll (1956). (Figure after Carroll) 
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Figure 15. Upper limits to fluorescent excitation rates assuming f and a 
equal 1. For the known f values of the red and violet systems, the curve 


is displaced accordingly. 








169 


Figure 16. (a) Coordinate system for treating the excitation of a polar molecule 
by electrons using the Born approximation (b) orientation of initial and final 


electron propogation vectors. 
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Figure 17. Temperature increase due to partical excitation by electrons and 
H I atoms. Each curve represents a constant observed excitation temperature 
so that the actual brightness is lower by this amount. Here we assume that H I - 
CN collisions predominantly satisfy electric dipole selection rules for CN exci- 


tation. 



HI DENSITY (cm- 3 ) 



ELECTRON DENSITY (cm' 3 ) 





Figure 18. of CN which would be expected foy excitation by a mixture of 

radiation and 125°K electrons or H I atoms as a function of observed T „ and 

10 

electron or H I density. Again we assume that the HI- CN cross section 
obeys dipole selection rules for CN excitation. 



HI DENSITY (cm -3 ) 



0 1.0 2.0 3.0 

ELECTRON DENSITY (cm' 3 ) 
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Figure 19. J = 0 -* 1 excitation cross-section for protons or singly charged 
ions as a function of velocity 


=1 


8k T 


7T m 


= 1. 31x 10 cm/sec 


P 


where T = 10 4 °K. 
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Figure 20. Coordinate system used for treating the excitation of a polar 
molecule by protons. 
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Figure 21 Excitation of the J = 1 -*1 and 0 ->2 transitions by protons as a function 

— 8 

of impact parameter. The proton velocity is v = 1. 342 v = 1. 75 x 10 cm/sec, 

where v is the most probable proton velocity 'at a kinetic temperature of 
4 

10 °K. The discrete points result from the numerical integration of Schrodinger's 


equation. 
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Figure 22. Excitation of the first-excited level of GH from the ground state by- 
protons as a function of impact parameter. As with Figure 21, the proton 
velcoity is 1. 75 x 10 6 cra/see. 



rr (cm sec 
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Impact Parameter, b (A ) 
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Figure 23. Number density of protons 'and electrons, as a function of the HII 

**■ 

kinetic temperature,, required to produce the 3°K excitation temperature of the 


J = 1 CN level. 
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Figure 24. Rotation temperature of the J = 2 and 3 CN levels, as a function 
of the H H kinetic temperature, when the number density of charged particles 
is such that the rotational temperature of the J = 1 level is 3°K. 
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Figure 25. Rotational temperature of the J = 1 CN level as a function of the 
number density of protons and electrons, for an HII temperature of 10, 000°K. 




N (cm" 3 ) 
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Figure 26. Velocity profiles of CN and Ca II absorption in the spectrum of 
Q Ophiuchi compared with H I 21cm emission in the same direction. 
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Figure 27 . Velocity profile of CN absorption m the spectra of X Persei and 
f Persei compared with H I 21cm emission in the corresponding directions. 
These two stars are only separated by about one degree of arc and appear to 


share the same interstellar cloud. 
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Figure 28. Velocity profiles of CN and Ca II absorption in the spectrum of 
20 Aquilae compared with H I 21cm emission in the same direction. 




HELIOCENTRIC VELOCITY (KM/SEC) 
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Figure 29: Comparison of the excitation temperature T^(CN) m Q Ophiuchi 

and C, Persei, with the series of direct measurements of the background rad- 
iation made by the Princeton group. 
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Figure 30. Three cases of error limit regions: 

A. Region is symmetric about most probable value which is 

< 0: area between B - rj and 8 + r? equals success 

probability. 

B. Region not symmetric about most probable value which is 

{ 0; area between zero and /J + 77 equals success probability 

C. Most probable value equals zero; area between zero and r\ 
equals success probability. 

In all cases, total area for positive /? is normalized to one. 







